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ABSTRACT 


The  Naval  Postgraduate  School  simulation  model,  QDYN86,  was  used  to 
examine  sputtering  of  nitrogen  from  the  (100)  faces  of  single  crystals 
of  molybdenum  and  tungsten.  The  nitrogen  placement  was  varied,  and 
analyses  were  conducted  on  the  sputtering  cross  sections  of  the 
nitrogen.  The  cases  where  the  adatom  was  directly  hit  by  the  incident 
ion,  or  if  it  was  sputtered  due  to  the  collision  cascade  process,  were 
analyzed  separately.  The  simulations  were  conducted  to  compare  the 
results  with  Winters’  recent  work,  and  to  build  upon  the  efforts  of  ear¬ 
lier  studies  completed  at  the  Naval  Postgraduate  School.  It  was  found 
that  placement  of  nitrogen  at  0.245  A  from  the  surface  of  molybdenum 
resulted  in  cross  sections  similar  to  those  found  by  Winters.  The 
effect  of  the  mass  of  the  substrate  was  verified,  in  that  a  substrate  of 
greater  mass  results  in  a  higher  sputtering  cross  section.  This  agreed 
with  Winters’  findings,  and  conflicted  with  earlier  conclusions  of  past 
theses.  The  adatoms  apparently  reduce  the  momentum  available  to  create 
collision  cascades,  reducing  the  sputter  yield  ratio  of  the  substrate 
when  the  ions  directly  hit  the  adatoms. 


3 


TABLE  OF  CONTENTS 


I .  BACKGROUND  7 

A.  HISTORICAL  OVERVIEW  7 

B.  PHYSICAL  UNDERSTANDING  10 

C  THEORETICAL  AND  EXPERIMENTAL  DEVELOPMENTS  15 

(1960  TO  PRESENT) 

D.  COMPUTER  SIMULATIONS  16 

E.  APPLICATIONS  18 

II.  OBJECTIVES  21 

A.  PREVIOUS  EFFORTS  25 

B.  WINTERS’  SINGLE  CRYSTAL  SPUTTERING  EXPERIMENTS  25 

C.  THESIS  OBJECTIVES  27 

III.  COMPUTER  SIMULATION  AND  MODEL  DEVELOPMENT  28 

A.  THE  COMPUTER  MODEL  AND  RELATED  PROGRAMS  28 

1  .  QDYN86  28 

2.  Ancillary  Programs  30 

B.  SUBSTRATE  AND  ADATOM  PROPERTIES  31 

C.  POTENTIAL  FUNCTIONS  33 

1.  General  33 

2.  Selection  of  Potential  Function  Parameters  35 

D.  TARGET  AND  IMPACT  AREA  CONCERNS  44 


4 


IV.  RESULTS  AND  DISCUSSION  47 

A.  GENERAL  47 

1.  Sputtering  Cross  Sections  47 

2.  Winters’  Results  47 

3.  Previous  Theses  at.  the  Naval  Postgraduate  School  48 

4.  Analyses  Conducted  in  This  Study  49 

B.  MOLYBDENUM  RESULTS  50 

t.  Nitrogen  Sputtering  Cross  Sections  and  Adatom  50 

P 1 acement 

2.  Nitrogen  Sputtering  Cross  Sections  as  a  55 

Function  of  Energy 

3.  Sputtering  of  the  Substrate  58 

C.  TUNGSTEN  RESULTS  58 

1.  Nitrogen  Sputtering  Cross  Sections  and  Adatom  58 

P 1 acement 

2.  Nitrogen  Sputtering  as  a  Function  of  Energy  66 

3.  Sputtering  of  the  Substrate  66 

D.  COMPARISON  OF  MOLYBDENUM  AND  TUNGSTEN  RESULTS  68 

1.  Comparison  of  Cross  Sections  68 

2.  Comparison  of  Substrate  Sputtering  Yields  70 

E.  DETERMINATION  OF  POSSIBLE  MASS  EFFECTS  70 

F.  COMPARISON  WITH  PREVIOUS  SIMULATIONS  74 

V.  CONCLUSIONS  AND  RECOMMENDATIONS  75 

APPENDIX  DETERMINATION  OF  SPUTTERING  CROSS  SECTIONS  77 

LIST  OF  REFERENCES  80 

INITIAL  DISTRIBUTION  LIST  85 


5 


ACKNOWLEDGEMENT 


I  wish  to  thank  those  who  preceded  me,  Steve  Webb  and  Dirk  Meyer- 
hoff,  and  to  acknowledge  the  fact  that  I  drew  quite  heavily  upon  their 
efforts.  1  also  wish  to  thank  Professor  Harrison  for  his  patience,  and 
the  concern  he  showed  throughout  this  exercise.  Special  thanks  to  my 
wife,  Lyn,  and  my  sons  Karl  and  Curtis  who  put  up  with  a  lot  during  the 
final  weeks  of  this  project.  Finally,  special  heartfelt  thanks  to  Dana 
Majors,  who’s  contribution  to  the  effort  could  not  go  unrecognized. 


6 


BACKGROUND 


A.  HISTORICAL  OVERVIEW 

When  a  surface  is  bombarded  by  a  beam  of  energetic  particles  under 
proper  conditions,  surface  damage  effects  can  be  observed.  This  damage 
can  manifest  itself  a  number  of  different  ways;  as  surface  damage  in  the 
form  of  pits,  blisters,  and  cones,  and  in  the  ejection  of  atoms  from  the 
target.  This  ejection  of  atoms  during  bombardment  is  called  sputtering. 

Sputtering  was  first  discovered  in  1853  by  Grove  [Ref.  1]  when  he 
observed  the  disintegration  of  the  cathodes  in  glow-discharge  tubes.  He 
noted  that  the  cathode  material  was  deposited  on  the  glass  walls  of  the 
tubes.  He  called  this  process  cathode  sputtering.  This  erosion  of  the 

cathode  was  also  noted  by  Plucker  [Refs.  2-4],  Gassiot  [Ref.  5]  and 
Faraday  also  reported  similar  observations. 

Fifty  years  passed  before  any  explanation  for  these  phenomena  was 
proposed.  In  1902  Goldstein  [Ref.  6]  presented  evidence  that  the  sput¬ 

tering  effects  reported  by  Grove  and  Faraday  were  caused  by  positive 
atoms  of  the  discharge  impacting  on  the  metal  cathode  of  the  tubes. 

Seven  years  later,  in  1909,  Stark  proposed  two  models  in  an  attempt  to 

explain  sputtering.  The  "hot  spot  model"  considered  the  sputtered  atoms 

to  be  the  result  of  evaporation  of  target  material  from  a  small  surface 
region  due  to  localized  heating  by  the  ion  beam.  The  "collision  model" 

proposed  that  sputtering  events  were  the  result  of  a  senes  of  binary 

collisions  initiated  by  a  single  ion  [Ref.  7].  These  models  served  to 

explain  the  experimental  results  at  the  time. 
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Later,  in  1921,  Thompson  [Ref.  6]  contributed  to  the  field  by 
suggesting  that  the  atomic  ejection  was  caused  by  the  release  of  radia¬ 
tion  as  the  ion  struck  the  target.  In  the  following  years,  a  consider¬ 
able  amount  of  theoretical  and  experimental  work  was  done.  Bush  and 
Smith  [Ref.  9]  in  1922  attempted  to  describe  the  ejection  of  atoms  as 

the  results  of  the  expansion  of  gas  adsorbed  by  the  target  material.  The 
following  year,  Kingdon  and  Langmuir  [Refs.  10,  11]  conducted  an  exper¬ 
iment  that  suggested  a  momentum  transfer  ejection  mechanism  for  sputter¬ 

ing.  They  bombarded  thoriated  tungsten  with  ions  in  a  glow  discharge 
tube.  This  was  a  special  case  of  sputtering,  where  the  thin  surface  film 

of  thorium  sputtered  rather  than  the  tungsten  substrate. 

In  1926,  von  Hippie  and  Blechschmidt  [Refs.  12-15]  proposed  a 
theory  that  described  sputtering  as  an  evaporation  of  the  surface  atoms. 
Through  spectroscopic  techniques,  von  Hippie  found  that  some  of  the 
sputtered  atoms  were  in  an  excited  state.  He  made  further  refinements 
on  Stark’s  hot  spot  model,  and  made  the  first  attempt  to  formulate  a 
sputtering  theory  on  the  basis  of  local  heating. 

Lamar  and  Compton  [Ref.  16],  in  1934,  published  "A  Special  Theory 
of  Cathode  Sputtering,"  which  led  to  the  "thermal  spike"  concept.  They 
suggested  that  binary  collision  processes  were  dominant  in  light -ion 

sputtering,  where  local  evaporation  predominates  for  heavy-ion  sputter¬ 
ing.  The  thermal  spike  concept  was  based  on  momentum  transfer  between 
the  bombarding  ion  and  the  lattice  atoms.  The  theory  suggested  that  a 
long  lived  high  temperature  volume  persisted  in  the  target  after  the 
collision  cascade  was  completed. 
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One  of  the  major  problems  facing  experimentalists  was  the  lack  of 


ability 

to 

reproduce  results. 

Penning  and  Moubis  in 

1940  published 

the 

results 

of 

their  studies  of 

the  effect  of  pressure 

on  sputtering 

yield 

[Ref.  17),  (the  sputtering  yield  is  defined  as  the  ratio  of  "sputtered" 


target 

atoms 

per 

incident  ion). 

They  found  that  the  sputtering 

yield 

was  reduced 

with 

an 

increase 

in  pressure. 

Collisions  between 

the 

ejected 

atoms 

and 

the 

surrounding 

gases  at  the 

higher  pressures  resulted 

in  the  backscattering  of  ejected  atoms  back  onto  the  target  surface, 
reducing  the  sputtering  yield.  When  the  background  pressure  was  kept 
below  10“5  Torr,  they  were  able  to  obtain  reproducible  results  for  ion 
energies  in  excess  of  500  eV. 

Keywell  made  an  attempt  to  formulate  Stark’s  collision  model  in 
terms  of  a  neutron  transport  model  originally  developed  for  nuclear  work 
[Ref.  18).  Wehner  [Refs.  19  -  21),  one  of  the  major  contributors  to  the 

modern  understanding  of  sputtering,  began  publishing  his  findings  in 
1954.  He  discounted  the  evaporation  model,  and  presented  strong  evi¬ 
dence  for  a  momentum  transfer  process.  He  demonstrated  the  effects  of 
crystal  structure  on  the  yield,  and  it  became  apparent  that  local  heat¬ 
ing  alone  could  not  account  for  the  effects  of  sputtering.  His  observa¬ 
tions  revived  interest  in  collision  theory.  One  of  his  major  contribu¬ 
tions  was  the  discovery  of  the  anisotropic  ejection  patterns  from  mono- 


crystalline  targets,  the 

now 

well 

known 

"Wehner  spot  pattern." 

In  1956  Harrison 

[Ref. 

22] 

applied 

statistical  methods  to 

sputter 

ing.  He  developed  a 

theory 

involving 

the  interaction  of  two 

distnbu 

tion  functions;  one  for  the  crystal  lattice  and  one  for  the  ion  beam. 
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He  utilized  transport  theory  ond  introduced  the  idea  of  ideal  collision 
cross  sections  to  explain  sputtering  from  amorphous  materials. 

Sputtering,  to  this  point,  had  been  regarded  as  an  annoying  mani¬ 
festation,  one  that  eroded  cathodes  and  filaments,  contaminated  plasmas, 
and  was  a  general  nuisance.  There  was  great  difficulty  in  obtaining 
reproducible  experimental  results.  The  theories  proposed  had  not  com¬ 
pletely  described  the  phenomena,  and  indeed  have  not  yet  today.  Before 
continuing  with  the  discussion  of  the  current  trends  in  sputtering 
theory  and  experimentation,  a  brief  discussion  of  the  currently 
’’accepted"  concepts  of  sputtering  is  in  order. 

B.  PHYSICAL  UNDERSTANDING 

Sputtering  is  the  ejection  of  atoms  from  a  target  when  bombarded  by 
energetic  ion  projectiles.  The  incoming  ion  collides  with  atoms  in  the 
bulk  of  the  target  material,  transferring  part  of  its  kinetic  energy 
and  momentum  to  the  target  atoms.  If  the  energy  transferred  to  an  atom 
is  greater  that  the  binding  energy  at  the  lattice  site,  then  a  pr imar  y 
r  eco i 1  a  tom  is  created.  This  recoil  a  tom  may  then  collide  with  other 
atoms  in  the  lattice,  creating  a  collision  cascade ,  transferring  energy 
throughout  the  material.  A  surface  atom  may  be  ejected  (sputtered)  if 
the  normal  component  of  the  energy  transferred  to  it  is  greater  than  the 
surface  binding  energy  of  the  lattice  [Ref.  23].  Figure  1  illustrates 
three  concepts  of  the  sputtering  process.  Figure  1  (a)  shows  the  inter¬ 
action  of  the  incident  ion  beam  with  the  surface  layer  of  the  target. 
Figure  1  (b)  illustrates  the  concept  of  the  "thermal  spike,"  where  a 
high  local  temperature  results  in  the  evaporation  of  surface  atoms. 
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cascade  process, 


which 


Finally,  Figure  1  (c)  illustrates  the  collision 

is  the  current  conceptual  model  for  sputtering  [Ref.  24]. 


Ion  Beam 


(a)  Sputtering  of  surface  atoms  as  a  result  of  the  inter¬ 
action  of  the  incident  ion  beam  and  the  target  sur¬ 
face. 

(b)  Surface  atoms  evaporated  as  a  result  of  a  thermal 
spike. 

(c)  Sputtering  from  a  collision  cascade,  where  energy  is 
directed  toward  the  surface  from  multiple  binary 
collisions.  [Ref.  24] 

Figure  1.  Models  of  the  Sputtering  Process. 
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The  mechanisms  for  the  sputtering  of  multicomponent  materials 


(targets  consisting  of  more  than  one  atomic  component)  are  more  compli¬ 
cated.  One  of  the  simpler  cases  is  the  situation  where  a  thin  layer  of 
a  third  component  is  placed  on  the  surface  of  the  thick  target  materia). 


This 

is  similar  to  the  case  examined  by  kingdon 

and 

Langmuir  in 

Refer- 

ence 

11. 

In 

this  case,  the  added  atoms 

on  the 

surface  are 

called  adatoms 

and 

the 

bulk 

of  the  material  is  called 

the  substrate 

[Ref. 

25]. 

Examples 

of 

this 

situation  were  covered  by 

Garrison, 

Wmograd, 

and 

Harrison 

[Ref. 

26] 

for 

oxygen  on  copper,  and  by 

Winters 

and 

Sigmund 

for 

nitrogen 

on  tungsten 

[Ref.  27]. 

Winters  proposed  three  mechanisms  for  the  sputtering  of  the  surface 
atoms,  illustrated  below  in  Figure  2.  Layer  (3)  is  the  thin  layer  of 

adatoms  on  the  surface  of  the  thick  target  (2).  An  ion  (1)  bombards  the 
target  surface.  In  Figure  2(a),  the  ion  hits  the  adatom,  and  the  adatom 
is  reflected  off  the  substrate,  either  directly,  or  after  it  penetrates 
the  target  slightly.  Figure  2(b)  illustrates  the  case  where  the  ion 

does  not  hit  the  adatom  directly,  but  penetrates  the  target,  and  is  in 

turn  reflected  itself.  The  reflected  ion  then  hits  the  adatom,  which  is 
then  sputtered  off.  Figure  2(c)  shows  the  case  where  the  ion  causes  an 
outward  collision  cascade,  much  like  that  illustrated  in  Figure  1(c). 

The  outward  flux  of  sputtered  substrate  in  turn  sputters  the  adatom. 

The  basic  processes  associated  with  sputtering  are  similar  to  those 


causing 

radiation 

damage  in 

the  bulk  of 

a  solid.  Sputtering, 

however, 

usually 

involves 

atoms 

in  the 

distorted  surface  layers  (selvage) 

of 

the 

target 

material. 

The 

term, 

"selvage,"  is 

derived  from  the  term 

"sel- 

vedge," 

which 

is  a 

narrow 

band  woven 

such  that  the  edge 

will 

not 
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unravel.  This  analogy  was  extended  to  the  distorted  layers  of  the  sub¬ 
strate  material  at  its  surface  [Ref.  28]. 


(a) 


(b) 


(c) 


Figure  2.  Three  Proposed  Mechanisms  for  the  Sputtering  of  a  Thin 
Layer  on  a  Thick  Target  [Ref.  27]. 

Sputtering  events  take  place  in  conditions  far  from  thermal  equili¬ 
brium,  and  are  not  evaporation  of  the  material.  Several  factors  can 
influence  the  sputtering  yield,  including  incident  ion  energy,  angle  of 
incidence,  ion  type,  target  material,  target  crystal  orientation,  the 
presence  of  adsorbed  molecules  on  the  surface  of  the  target,  ambient 
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pressure,  and  numerous  other  factors. 


In  order  to  obtain  reproducible 


experiments,  the  following  conditions  must  be  met  [Ref.  29]: 

1.  The  target  surface  must  be  clean,  that  is  free  of  contaminants  in 
the  form  of  adsorbed  gases,  lubricants  from  the  vacuum  pump, 
cleaning  solvents,  etc. 

2.  The  gas  pressure  must  be  such  that  the  mean  free  path  of  ions  and 
sputtered  atoms  is  large. 

3.  The  ion  current  density  must  be  high  and  the  background  pressure 

low  so  that  formation  of  surface  layers  is  prevented  during  the 
experiment. 

4.  The  ions  must  strike  the  target  at  a  known  angle. 

5.  The  energy  spread  of  the  incident  beam  must  be  small. 

6.  The  ionizing  conditions  in  the  ion  source  should  be  such  as  to 

minimize  the  production  of  multiply  charged  species;  the  ions 

must  be  uniformly  charged  and  mass  separated. 

7.  The  lattice  orientation  of  monocrystalline  targets  must  be  known. 

The  field  of  sputtering  is  rich  with  terms  used  to  describe  the 

various  observed  and  theoretical  aspects  of  sputtering.  Specific  terms 

will  be  introduced  as  necessary  to  describe  certain  events,  but  two 

terms  do  merit  mention  here.  T  ransmission  sputtering  is  the  ejection  of 
atoms  from  the  rear  of  a  thin  target.  This  occurs  when  sufficient 

energy  is  transported  throughout  the  target  to  allow  atoms  to  overcome 

the  potential  energy  binding  it  to  the  target,  allowing  it.  to  escape. 

Back-sputtering  is  the  more  familiar  form,  where  atoms  are  ejected  from 

the  surface  of  the  target  material.  Two  other  terms  widely  used,  espe¬ 
cially  when  discussing  the  sputtering  of  multicomponent  materials  are 

physical  sputtering  and  chemical  sputtering.  Physical  sputtering 
involves  a  transfer  of  kinetic  energy  from  the  incident  particle  to  the 
atoms  in  the  target,  and  the  subsequent  ejection  of  the  atoms.  Back 
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sputtering  and  transmission  sputtering  are  manifestations 


o  f  physical 


sputtering.  Chemical  sputtering  is  the  result  of  a  chemical  reaction 
induced  by  the  bombarding  particle,  which  produces  an  unstable  chemical 
compound  on  the  surface  of  the  target  [Ref.  30]. 

C.  THEORETICAL  AND  EXPERIMENTAL  DEVELOPMENTS  (1960  TO  PRESENT) 

In  1962,  Wehner  and  others  [Ref.  31]  published  data  on  the  sput¬ 
tering  yields  of  metals  and  semiconductors  in  the  100-600  eV  range. 
Later  in  1966,  Wehner  [Ref.  32],  in  a  report  for  Litton  Systems  pub¬ 
lished  detailed  results  covering  years  of  low  energy  sputtering 
research. 


Silsbee  [Ref. 

33],  in 

an 

effort 

to 

account  for  the 

angular 

distri- 

bution  of 

ejecta, 

the  "Wehner 

spot 

patterns,"  proposed  a 

focused 

colli- 

sion  model, 

which 

allowed 

the 

transport 

of  momentum  in 

cry s  tals 

along 

preferred 

directions. 

Experimental 

results  had  indicated 

that 

the 

sputtering 

yield 

for 

single  crystal 

targets  were  dependent 

upon 

the 

crystallographic 

orientation  of  the 

target  and  the  incident 

ion 

beam. 

This  seemed  reasonable,  considering  the  "holes"  observed  in  crystal 
models  of  when  viewed  from  different  planes.  Focusing  can  be  seen  to 
contribute  to  the  development  of  collision  cascades  within  the  target, 
it  could  not  explain  the  Wehner  spot  patterns. 

Sigmund  and  Lehmann  [Ref.  34]  proposed  an  alternative  mode),  based 
on  Boltzman  transport  equations,  requiring  the  target  surface  to  have  an 
ordered  surface.  In  this  model,  an  atom  in  the  collision  cascade 
would  sputter  if  its  kinetic  energy  component  normal  to  the  surface  was 


15 


sufficient  to  overcome  the  surface  potential  barrier.  Thompson  [Ref. 


35]  proposed  another 

model  to 

account 

for 

the 

spot  patterns, 

in  which 

the  surface  attraction 

for 

the 

escaping 

atom 

causes  a  refraction  of  its 

velocity  vector  away 

from 

the 

normal, 

resulting 

in  a  distortion 

of  the 

angular  distribution  of 

ejecta. 

The  theories 

of 

Sigmund  and 

Thompson 

predict  the  sputtering  yields  of  amorphous  or  polycrystalline  targets 
fairly  well,  but  do  not  accurately  reflect  the  experimental  results  from 
the  ordered  surfaces  of  single  crystal  targets. 

Numerous  theories  have  been  proposed,  but  none  succeeds  in  fully 
describe  all  aspects  of  sputtering  phenomena.  Harrison  [Ref.  36]  in  a 
recent  review,  stated  that  there  were  currently  no  less  than  seven  types 
of  sputtering  theories  in  the  literature.  The  statistical  theories  of 
Thompson  [Ref.  35]  and  Sigmund  [Ref.  34]  predict  the  yield  relatively 

well,  and  provide  information  on  the  ejected  atom  energy  distribution 
function,  and  ejected  atom  angular  distribution  for  polycry  stallme 
targets.  Kelly  [Ref.  37]  classifies  sputtering  processes  according  to 

time  scales. 

During  this  time,  computer  simulations  were  used  to  study  sputter¬ 
ing.  The  models  and  codes  now  in  use  have  been  evolving  for  the  past  30 

years.  The  next  section  will  present  some  basic  philosophy  of  computer 
simulations,  with  applications  to  sputtering  research, 

D.  COMPUTER  SIMULATIONS 

In  1960  a  new  tool  was  added  to  the  scientific  arsenal,  aimed  at 

examining  sputtering.  Gibson,  Goland,  Milgram  and  Vineyard  [Ref.  36] 
published  a  report  of  the  first  computer  simulation  examining  radiation 


16 


damage  in  a  material.  They  simulated  metallic  copper  and  studied  radia¬ 
tion  damage  events  at  low  and  moderate  energies,  up  to  400  eV.  In  the 

model,  one  atom  was  given  an  arbitrary  Kinetic  energy  and  direction  of 

motion,  simulating  its  having  been  struck  by  an  energetic  particle. 
This  was  one  of  the  first  published  accounts  of  a  computer  simulation  to 

study  such  events. 

The  high  operating  speed  of  computers  makes  them  the  natural  choice 

for  processing  the  numerous  calculations  required  in  a  numerical  analy¬ 
sis.  The  use  of  a  computer  simulation  frees  one  from  the  constraints  of 

general  theories,  and  allows  one  to  examine  the  basic  physics  of  the 

system,  to  see  how  theory  and  experiment  interact.  In  the  words  of 

Harrison  [Ref.  36,  p.  4]: 

A  simulation  is  not  a  theory:  it  is  a  mathematical  tool  which  is  used 
to  test  the  fundamentals  of  a  theory.  The  computer  can  model  a 

system  with  a  minimum  set  of  physical  assumptions.  This  inherent 

simplicity  helps  to  elucidate  complex  problems  like  sputtering: 
Simulations  develop  ideas  which  can  then  be  exploited  by  both 

experimentalists  and  theorists. 

Computer  simulations  fall  into  two  general  categories:  time-step 

models  and  event-store  models.  An  event-store  program  moves  from  event 
to  event,  skipping  the  intervening  time.  It  maintains  a  list  of  poten¬ 
tial  future  events  (hence  the  name  "event-store"),  and  checks  and 
updates  the  lists  to  determine  which  event  happens  next.  This  model 
works  well  when  the  model  is  well  understood,  and  the  events  are  well 
separated  in  time.  A  time-step  program  carries  the  model  forward  for  a 
short  period  of  time,  computes  everything  that  happens  to  the  system  in 
that  time  interval,  updates,  and  then  continues  on.  This  program  is 

most  useful  when  several  things  happen  simultaneously.  The  time-step 
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program  tends  to  be  shorter  than  the  event  store  program,  but  it  also 


runs  slower. 

The  event-store  model  used  in  sputtering  research  is  based  on  the 
binary  collision  approximation  (BC)  [Ref.  39].  The  assumption  made  is 
that  each  particle  only  interacts  with  one  other  particle  at  a  time,  and 
this  other  particle  is  usually  assumed  to  be  stationary.  These  models 
are  inherently  linear  calculations  [Ref.  40]. 

The  time-step  model  used  in  this  investigation  is  based  on  simulta¬ 
neous  multiple  interactions  (Ml).  Newtons  laws,  usually  expressed  in 
Hamiltonian  form,  are  numerically  solved  for  many  particles.  The  Ml 
model  used  in  this  thesis,  QDYN86,  is  the  latest  revision  of  the  Ml 
program  developed  and  used  at  the  Naval  Postgraduate  School,  and  Penn¬ 
sylvania  State  University.  QDYN86  is  a  full-lattice  simulation,  which 
models  the  dissipation  of  an  incident  ion’s  momentum  in  a  single  crystal 
target,  using  classical  mechanics.  This  program  can  generate  different 
surfaces  of  several  crystal  structures.  Adatoms  can  be  placed  on  t.he 
surface  of  the  crystal  to  simulate  reacted  surfaces.  Specific  details  of 
the  simulation  will  be  presented  m  a  later  section  of  this  thesis. 

As  can  be  seen,  a  considerable  effort  has  been  expended  to  study 
the  mechanisms  of  sputtering.  There  is  a  practical  reason  for  the 
interest,  beyond  the  lure  of  pure  research. 

E.  APPLICATIONS 

Sputtering  was  long  regarded  as  an  undesired  and  little-understood 
effect  [Refs.  41,  42].  It  destroyed  cathodes  and  grids  in  gas  discharge 
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tubes,  and  contaminated  plasmas  and  the  surrounding  walls  [Ref.  43], 


There  was  great  concern  about  damage  to  spacecraft  and  satellites  from 
sputtering.  Efforts  to  understand  and  reduce  these  effects  provided 
major  impetus  to  the  study  of  sputtering. 

However;  sputtering,  in  the  form  of  the  controlled  removal  of 


surface 

atoms  from  a  target,  is 

becoming 

especially  important 

in 

its 

own 

right. 

The  ability  to  precisely 

control 

an  ion  beam,  and 

to 

remove 

atoms 

from  surfaces  is  a  very 

important 

tool  in  research, 

and 

in 

the 

manufacture  of  miniature  components. 

Sputter  ion  sources  can  be  used  for  cleaning  of  surfaces,  by  remov¬ 
ing  adsorbed  surface  molecules  to  a  degree  that  is  impossible  to  achieve 
chemically  or  mechanically.  The  atoms  sputtered  from  a  surface  can  be 

analyzed  in  a  mass  spectrometer,  giving  information  about  the  surface 
composition.  One  of  the  more  important  commercial  applications  of 
sputtering  is  the  deposition  of  thin  films  on  a  large  variety  of  sub¬ 
strates,  especially  useful  in  the  manufacture  of  microelectronics. 

Sputtering  is  used  in  micromachining  and  depth  profiling  of  thin  films 
[Ref.  23]. 

There  are  applications  of  erosionally  modified  surfaces  in  the 
fabrication  of  optical,  magnetic  and  surface  acoustical  technologies. 
Areas  include  grating  fabrication,  magnetic  bubble  technology,  ion 
polishing,  and  reactive  ion-beam  etching  [Ref.  44],  Research  is  being 
conducted  to  examine  possible  biomedical  applications  of  sputtering. 
These  include  surface  modification  of  biomedical  materials,  using  sput¬ 
tering  for  pathological  discrimination,  and  with  applications  to 

implants  and  prostheses  [Ref.  45]. 
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One  of  the  first  commercial  applications  of  sputtering  was  the  use 


of  sputtering 

for 

the  deposition 

of 

solid  film  lubricants  [Ref. 

46]. 

Sputtering  had 

the 

advantage  that 

it 

allowed  for  the  deposition 

of  a 

variety  of  materials,  on  a  large  variety  of  substrates.  The  sputtered 
films  are  very  dense,  and  strongly  adherent.  These  properties  are 
particularly  useful  for  corrosion  resistance  and  lubrication. 

The  list  of  possible  applications  could  continue.  The  ability  to 
understand  and  control  the  surface  characteristics  of  a  material  to  the 
atomic  level  is  a  very  powerful  tool.  The  majority  of  the  work,  espe¬ 
cially  m  regards  to  the  applications,  is  the  result  of  considerable 
experimental  effort.  If  one  could  better  understand  the  mechanisms 
underlying  the  observed  manif  estations,  then  experiments  could  be  better 
designed,  and  applications  could  be  more  quickly  realized.  It  is  toward 
this  end  that  this  research  is  conducted. 
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II. 


OBJECTIVES 


A.  PREVIOUS  EFFORTS 

The  interaction  of  nitrogen  with  the  (100)  plane  of  tungsten  has 
been  widely  examined.  It  was  considered  to  be  a  useful  system  on  which 
to  base  the  development  of  adsorption  kinetic  and  dynamic  models 

[Ref.  47].  Winters  [Ref.  48],  m  1982,  published  a  paper  comparing  the 
sputtering  of  chemisorbed  nitrogen  from  poly  crystal  line  targets  of 
molybdenum  and  tungsten.  His  experiments  showed  that  the  nitrogen 
sputtering  yield  tended  to  increase  as  the  atomic  weight  of  the  target 

substrate  increased. 

This  paper  of  Winters  provided  the  basis  for  two  Masters  theses  at 

the  Naval  Postgraduate  School.  In  1983,  Meyerhoff  [Ref.  49]  used  an 
earlier  version  of  Harrison’s  computer  simulation  sputtering  model  to 


study  sputtering  from 

nitrogen 

reacted 

(100)  surfaces 

of  tungsten  and 

molybdenum  targets. 

His  goal 

was  to 

compare 

the 

simulation  results 

with  the  experimental 

results 

recently 

published 

by 

Winters  in  Refer- 

ence  48.  Meyerhoff  concentrated  on  the  reported  mass  effects.  He 
concluded  that  while  there  may  be  a  mass  relation,  the  distance  of  the 
adatom  from  the  substrate  had  a  much  more  profound  effect  on  the  sput¬ 
tering  yield. 

In  1986,  Webb  [Ref.  50]  sought  to  further  examine  the  adatom  pla¬ 
cement  problem.  Webb  again  used  (100)  surfaces  of  monocrystalline 
tungsten  and  molybdenum  reacted  with  nitrogen  and  bombarded  with  argon 
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ions.  He  narrowed  the  range  of  the  placement  of  the  adatoms,  and  his 


results  indicated  that  the  adatoms  were  positioned  above  the  (100)  sur¬ 
face  of  the  target  plane. 

In  all  these  cases,  the  nitrogen  adatoms  were  assumed  to  be  in  the 
four-fold  position,  as  shown  below  in  Figure  3.  heyerhoff  assumed  the 
nitrogen  atom  was  slightly  above  the  surface  of  the  target  plane,  in  a 
position  of  equal  distance  from  its  five  nearest  neighbors.  Webb  exa¬ 
mined  the  nitrogen  sputtering  yield  with  the  nitrogen  placed  at  three 
different  heights  for  molybdenum  and  two  heights  for  tungsten. 


Figure  3.  Nitrogen  Adatom  Placed  in  Four-fold  Position  in  BCC 

Unit  Cell. 

The  basis  for  the  four-fold  placement  assumption  is  well  estab¬ 
lished  in  the  literature.  An  earlier  work  by  Clavenna  and  Schmidt 

examined  the  interaction  of  Ng  with  W(100)  [Ref.  51).  They  studied  the 
binding  states  and  condensation  and  desorption  kinetics  of  nitrogen  on 

W(100)  using  flash  desorption  spectrometry.  They  felt  that  the 
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nitrogen  was  in  the  four -fold  position,  but  were  uncertain  exactly 


where  the  nitrogen  was  placed  vertically  in  relation  to  the  target  sur¬ 
face.  Adams  and  Germer  [Refs.  52,  53]  also  examined  the  adsorption  of 

nitrogen  on  tungsten.  They  used  a  combination  of  Low  Energy  Electron 
Diffraction  (LEED),  f  lash-desorption  mass-spectrometry  and  contact 
potential  measurements  in  their  experiments.  They  refined  their  exper¬ 
iment  by  using  the  15N-isotope  of  nitrogen,  in  order  to  distinguish 

between  nitrogen  and  CO  present  in  the  background  gas.  They  felt  that 
since  the  nitrogen  atom  was  half  the  size  of  the  tungsten  atom,  that  it 
would  sit  in  the  well  formed  by  the  surface  atoms  of  the  substrate. 

Griffiths,  Kendon,  King  and  Pendry  [Ref.  54]  examined  target  atom 
displacement  due  to  the  introduction  of  the  nitrogen  adatoms.  They 
performed  as  series  of  LEED  experiments  with  various  nitrogen  coverages 
and  beam  energies.  They  found  that  a  fractional  coverage  of  approxi¬ 

mately  0.4  (defined  as  the  ratio  of  the  adatoms  to  the  number  of  sub¬ 
strate  atoms  for  a  specific  face  of  the  target)  was  relatively  stable. 
At  higher  coverages,  of  about  0.5,  the  nitrogen  would  be  absorbed  into 
the  bulk  at  temperatures  of  about  1000  K.  The  0.4  coverages  were 

stable  at  this  temperature.  The  model  suggested  by  Griffiths,  et.  al. 
was  a  contracted  domain  structure,  where  the  nitrogen  occupies  a  four¬ 

fold  hollow  site,  but  that  the  four  surface  tungsten  atoms  are  uni¬ 
formly  displaced  towards  the  adatom.  This  leads  to  a  series  of  "is¬ 

lands"  consisting  of  16  nitrogen  atoms  and  the  corresponding  surface 
tungsten  atoms,  as  shown  in  Figure  4.  This  type  of  formation  was  phys¬ 
ically  modeled  by  conducting  a  series  of  laser  diffraction  experiments. 
Laser  diffraction  gratings  were  constructed,  with  an  average  island 
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size  of  4x4,  and  when  illuminated  by  laser  light  produced  diffraction 
patterns  similar  to  those  of  the  LEED  beams  on  the  actual  target. 


Again,  while 

this 

provided  insight  into  the 

relative  placement 

of 

the 

nitrogen  on 

the 

surface,  it  did  not  provide 

any  information  on 

how 

the 

nitrogen  was  related  vertically  to  the  crystal  surface. 


Illustrates  contracted-domain  structure  of  the  0.4 
monolayer  of  Nitrogen  on  W(100)  or  Mo(IOO)  surface. 
Large  hatched  circles  illustrate  the  top  layer  of  the 
substrate  atoms,  the  small  filled  circles,  the 
nitrogen  atoms. 


a) 

Plan  view:  Shows  Domain 

and 

Boundary  of 

Structure. 

b) 

Cross 

Section  through 

line 

AA.  [Ref. 

54,  p.1586]. 

Figure  4.  Suggested 

"Island"  Formation 

of 

Nitrogen  on 

T  ungsten. 
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Meyerhoff  assumed  the  equilbnum  position  of  the  nitrogen  in  his 


research.  Webb  examined  equilibrium  points  both  above  and  below  the 

substrate  surface  plane.  While  there  was  a  great  body  of  information 
available  on  the  adsorption  of  nitrogen  on  single  crystal  surfaces, 

there  was  little  data  on  the  sputtering  of  nitrogen  from  single  crystal 
surfaces.  This  was  rectified  by  a  recent  report  by  Winters  [Ref.  55]. 

B.  WINTERS’  SINGLE  CRYSTAL  SPUTTERING  EXPERIMENTS 

This  paper  [Ref.  55]  was  a  continuation  of  Winters’  previous 

studies  in  the  sputtering  of  chemisorbed  nitrogen  on  tungsten.  In  his 
1974  paper  he  proposed  two  mechanisms  for  the  sputtering  of  nitrogen 

from  the  tungsten  surface  [Ref.  27].  He  suggested  that  for  low  energy 
incident  ions,  the  primary  mechanism  for  the  sputtering  of  the  nitrogen 
was  direct  knock-on  collisions  with  the  impinging  or  reflected  bombard¬ 
ing  ions.  The  proposed  mechanism  for  sputtering  at  higher  bombarding 
energies  included  nitrogen  atoms  knocked  away  by  sputtered  substrate 

atoms.  Winters’  1982  paper  [Ref.  48]  examined  the  sputtering  of 
nitrogen  from  molybdenum  and  tungsten  polycrystalline  targets. 

This  latest  report  by  Winters  and  Taglauer  described  experiments 
investigating  the  sputtering  of  chemisorbed  nitrogen  from  W(100), 
W(  1 11),  W(110)  and  Mo(100)  single  crystal  surfaces.  The  targets  were 

bombarded  with  helium,  argon  and  xenon  ions  in  the  energy  range  of  300 
eV  to  5000  eV.  The  relation  between  the  s  or  bate  mass,  substrate  mass 
and  sputter  yield  were  examined.  Conclusions  were  drawn  about  thermal 
spikes,  recoil  implantation,  and  cascade  mixing.  Since  these  exper¬ 
iments  were  conducted  on  single  crystal  surfaces  they  provided  some 
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basic  information  on  the  physics  of  the  sputtering  of  multicomponent 


systems.  Additionally,  since  the  experimental  data  was  derived  from 
single  crystals,  they  lend  themselves  more  directly  to  computer  model¬ 
ing.  A  brief  description  of  the  experiment  indicates  how  the  system 

can  be  modeled  in  the  computer  simulation. 

A  Faraday  cup  and  the  four  single  crystal  targets,  W(lOO),  W(111), 

W(110),  and  Mo(IOO)  were  mounted  on  a  rotating  manipulator.  The  tar¬ 

gets  were  cleaned  by  heating  to  about  2500K  for  tungsten  and  2200K  for 
molybdenum  by  RF  induction  heating  [Ref.  56].  The  sample  was  exposed 
to  until  about  1/2  monolayer  was  adsorbed.  This  is  in  agreement 

with  the  findings  of  stable  nitrogen  adsorption  on  tungsten  as  men¬ 
tioned  in  references  52  thru  54.  The  use  of  the  isotope  allowed 

one  to  isolate  the  adsorbed  species  involved  in  the  study. 

After  the  adsorption,  an  Auger  spectrum  was  taken  to  ensure  the 

absence  of  impurities  from  the  sample  surface.  The  ion  beam  was 
scanned  by  the  Faraday  cup  in  order  to  ensure  beam  uniformity,  and  to 

measure  the  beam  current  density.  The  ion  beam  was  rastered  across  a 
0.48  cm  aperture,  impinging  at  normal  incidence  upon  the  sample  cen¬ 
tered  behind  the  aperture.  A  nitrogen  Auger  peak-to-peak  intensity  was 

determined  as  a  function  of  the  ion  dose.  Finally  the  sample,  which 

was  partially  cleaned,  was  exposed  to  the  ambient  for  approximately  600 
sec  after  each  run  to  monitor  the  readsorption  of  nitrogen.  If  read- 

sorption  was  observed,  the  run  was  discarded. 

The  care  taken  in  these  experiments  is  quite  apparent.  The 
parameters  of  this  experiment  lend  themselves  quite  well  to  computer 

simulation.  Winters’  also  presented  a  comparison  of  the  experimental 
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results  with  theoretical  calculations  based  on  the  Sigmund-Winters 


theory  for  the  sputtering  of  chemisorbed  gas  [Ref.  27].  One  recommend¬ 
ation  in  the  paper  was  : 

In  particular,  it  would  be  useful  to  compare  an  extensive  set  of 
molecular  dynamics  calculations  with  both  the  experimental  data  and 
also  the  calculations  presented  in  this  paper.  [Ref.  55,  p.  26] 


C.  THESIS  OBJECTIVES 

The  primary  objective  of  this  thesis  is  to  examine  sputtering  from 
the  (100)  surfaces  of  single  crystals  of  molybdenum  and  tungsten. 
Nitrogen  will  be  placed  in  different  equilibrium  positions  in  relation 


to  the 

substrate  surface,  and  the 

system  will 

be 

bombarded 

with 

argon 

ions  of 

various  energies  at  normal 

incidence. 

The 

effect  of 

the 

mass 

of  the  substrate  in  relation  to  the  sputtering  yield  will  be  examined. 
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III.  COMPUTER  SIMULATION  AND  MODEL  DEVELOPMENT 


A.  THE  COMPUTER  MODEL  AND  RELATED  PROGRAMS 
1 .  QDYN86 

The  computer  simulation  used  in  this  study  is  called  QDYN86, 
which  stands  for  the  1986  revision  of  the  QDYN  (Quick  DYNamics)  program 
used  by  Harrison  at  the  Naval  Postgraduate  School.  This  program  uses 
multiple-interaction  (Ml)  logic  in  a  t ime-s tep  approach.  The  time-step 
logic  is  appropriate  when  many  events  occur  simultaneously.  The  model 
is  based  on  classical  mechanics,  using  Newton’s  laws  expressed  in  Hamil¬ 
tonian  form  to  simplify  the  calculations.  The  initial  inputs  to  the 
program  include  the  target  crystalline  structure,  atomic  masses  and 
potential  functions  of  the  substrate,  incident  ion,  and  any  adatoms; 
adatom  locations,  and  bombarding  ion  angle  of  incidence,  energy  and 
impact  point.  The  program  develops  the  subsequent  collision  cascade, 
tracking  the  positions  and  velocities  of  the  target  atoms  through  time. 
The  computations  terminate  when  there  is  insufficient  energy  for  any 
further  ejections  to  occur. 

The  positions  and  velocities  of  the  moving  atoms  are  tracked 
through  each  time-step.  In  order  to  strike  a  compromise  between  exces¬ 
sive  computer  run  time,  and  yet  to  maintain  reasonable  energy  conserva¬ 
tion,  the  time-step  increment,  is  variable.  The  time-step  is  determined 
by  a  specified  distance  divided  by  the  highest  atomic  velocity  in  the 
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target. 


The  distance  chosen  is  0.1  lattice  units  (LU) ,  where  the  lat¬ 


tice  umt  is  defined  as  one-half  the  lattice  parameter,  a0  (for  a  cubic 
structure).  The  time-step  increment  is  further  modified  by  other  fac¬ 
tors  which  take  into  consideration  the  previous  velocity,  smoothing  the 
transition  m  time. 

The  velocities  and  positions  of  every  atom  m  the  target  are 
not  calculated  for  each  atom  at  every  time-step.  Again,  this  would 
consume  excessive  amounts  of  computer  time.  As  a  result,  the  forces  are 
not  computed  on  a  specific  target  atom  until  the  atom  is  struck  by  a 
moving  atom.  Atoms  which  rise  above  the  target  surface  are  put  in  a 
tentative  list  of  ejected  atoms.  After  further  atom  ejection  from  the 
target  is  unlikely,  the  atoms  in  this  list,  are  tested  to  see  if  they 
were  travelling  with  sufficient  velocity  to  overcome  the  attractive 
forces  of  the  target  atoms.  Atoms  which  fail  this  test  are  added  to  the 
target,  and  are  not  counted  as  having  been  sputtered.  Specific  actions 
are  taken  for  atoms  that  leave  the  bottom  of  the  target,  and  through  the 
sides.  The  program  ma  i nta i ns ,  and  periodically  updates  a  listing  of  an 
atoms  nearest  neighbors.  After  the  completion  of  a  collision  cascade, 
the  program  reinitializes,  and  prepares  for  another  trajectory.  A  num¬ 
ber  of  trajectories  are  run  at  different  impact  points  in  order  to 
obtain  better  statistical  results. 

Subsequent  sections  in  this  chapter  will  provide  further 
details  on  the  choice  of  potential  functions,  substrate  characteristics, 
impact  points,  and  other  factors  as  they  pertain  to  this  study.  A  num¬ 
ber  of  other  articles  provide  more  detailed  information  on  the  QDYN 
simulation,  and  other  molecular  dynamics  simulations.  Reference  36 
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provides  a  good  overview  of  sputtering  models.  Harrison  and  Jakas  com¬ 


pared  the  HI  and  binary  collision  models  [Ref.  57],  and  with  Webb, 
examined  a  hybrid  code  between  the  two  models  [Ref.  56].  The  most 
recent,  and  most  detailed  discussion  of  the  QDYN  program  by  Harrison  and 
Jakas  identified  uses  of  the  program  for  studies  involving  insulators 
and  semiconductors  [Ref.  59]. 

2 .  Anc i I  I ar y  Programs 

A  number  of  other  programs  are  used  to  support  the  main  simu¬ 
lation.  These  programs  are  identified  and  are  briefly  described  below. 

a.  QDYNLIBB  (QD66L I BB) 

This  program  forms  a  library  from  which  the  mam  program 
derives  routines  for  calculating  the  forces  and  potentials.  The  poten¬ 
tial  and  forces  tables  are  dimensioned  to  10000,  and  support  four  poten¬ 
tial  functions.  Adatom  routines  are  available  for  placing  adatoms  on 
the  surface  of  the  target,  or  creating  steps  or  vacancies. 

b.  TARGLIBB 

This  program  forms  a  library  of  routines  used  for  generat¬ 
ing  (100),  (110),  and  (111)  faces  of  BCC,  FCC  and  Diamond  lattice  tar¬ 
gets.  A  rotated  FCC  (001)  lattice  can  also  be  generated. 

c.  ANMOL 

This  program  analyzes  the  basic,  data  set  generated  by  the 
ma  i  n  program.  The  output  of  the  main  program  consists  mainly  of  a  list¬ 
ing  of  the  positions  and  velocities  of  the  atoms.  ANMOL  uses  this  data 
file  to  calculate  which  atoms  are  sputtered,  from  which  layer  m  the 
target  lattice  they  originated,  and  rotational  and  vibrational  analysts 
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of  the  ejected  atoms.  Determination  of  multimers,  ejection  time  dis 


tr i but  ions  and  ejection  energies  are  also  calculated. 

d.  ANPLOT 

This  is  a  program  used  to  interface  with  the  DISSPLA 
graphics  package.  ANPLOT  is  used  to  generate  graphical  analyses  of  the 
data  file  generated  by  the  main  program. 

e.  POTTER 

This  program  is  used  to  generate  potential  functions.  An 
abbreviated  input  deck  is  read  by  POTTER,  and  two  data  files  are  gener¬ 
ated.  These  files,  FORCE  DATA  and  ENERGY  DATA,  are  used  by  the  next 
program  to  obtain  a  graphical  representation  of  the  functions. 

f.  POTPLOT 

This  is  another  graphics  routine  that  interfaces  with  the 
DISSPLA  graphics  package.  POTPLOT  uses  the  data  files  generated  by 
POTTER,  and  produces  a  liner  plot  of  the  potentials,  and  a  semi-log  plot 
of  the  forces.  This  allows  one  to  quickly  determine  which  modifications 
must  be  made  to  obtain  the  proper  shaped  curve  for  the  desired  potential 
and  force  function.  The  potential  functions  used  in  this  thesis  were 
generated  using  this  program. 

B.  SUBSTRATE  AND  ADATOM  PROPERTIES 

The  systems  of  nitrogen-on-molybdenum  and  nitrogen-on-tungsten  were 
chosen  by  Winters  for  his  sputtering  experiments  for  a  very  specific 
reason.  Many  of  the  physical  parameters  of  the  two  substrate  materials 
are  very  similar,  and  the  behavior  of  nitrogen  with  each  was  thought  to 
be  almost  identical  also.  This  provided  a  means  of  physically  reducing 
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the  variables  m  the  system,  enabling  him  to  isolate  and  concentrate  on 


the  mass  difference  of  the  two  materials.  A  surrmary  of  the  Key  physical 
parameters  for  the  two  systems  is  show  below  in  Table  1. 

TABLE  1.  PHYSICAL  DATA  FOR  MOLYBDENUM  AND  TUNGSTEN. 

MOLYBDENUM  TUNGSTEN 


Atomic  Number  (Z) 

Atomic  Weight  (amu) 
Atomic  Radius  (A) 

Density  (gm/cm^) 

Crystal  Type 
Lattice  Constant  (A) 
Lattice  Unit  (1/2  a0,  A) 
Re  (A) 

Va I ence 

ionic  Radius  (A) 

Cohesive  Energy  (eV) 
Binding  Energy  of  N  (eV) 


42 

74 

95.94 

183.85 

1 . 363 

1 .371 

10.22 

19.3 

BCC 

BCC 

3. 147 

3.  165 

1 .5735 

1 . 5825 

2.8 

2.894 

+  4 

+4 

0. 70 

0.70 

6.82 

8.90 

6.5 

6.5 

Data  for  Table  1  derived  from  References  60-62. 


Table  1  illustrates  the  high  degree  of  similarity  between  some  of 
the  physical  properties  of  molybdenum  and  tungsten.  Nitrogen  is  assumed 
to  occupy  the  four-fold  position  on  the  surface  of  both  materials.  The 
major  difference  between  the  two  materials  is  the  significant  difference 
in  mass  and  Z.  The  pertinent  physical  properties  of  the  adatom, 
nitrogen;  and  the  incident  ion,  argon;  are  summarized  below  in  Table  2. 

TABLE  2.  PHYSICAL  DATA  FOR  NITROGEN  AND  ARGON. 

N i trogen  Argon 

Atomic  Number  (Z)  7  18 

Atomic  Weight  (amu)  15  39.94 

Data  for  Table  2  derived  from  Reference  63. 
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c. 


POTENTIAL  FUNCTIONS 


1 .  Genera  I 

The  dynamics  of  t.he  atomic  interactions  are  controlled  by  the 
choice  of  the  potential  functions  used  to  describe  the  potential 
energies  and  forces  felt  by  the  atoms.  There  is  no  "right  choice"  for  a 
potential  function  in  these  types  of  simulations,  indeed,  Harrison  has 
called  the  development  of  useful  potential  functions  for  simulations  a 
"black  art".  A  number  of  different  potential  functions  can  and  have 
been  used  in  sputtering  simulations.  Harrison  has  found  that  the  choice 
of  potential  function  parameters  do  not  overly  affect  the  outcome  of  the 
simulation  [Ref.  64).  The  potential  functions  used  in  this  investiga¬ 
tion  are  described  briefly  below.  Detailed  discussions  of  these,  and 
other  potential  functions  can  be  found  in  Torrens  [Ref.  65).  A  number 
of  survey  articles  by  Harrison  describe  more  fully  how  the  potential 
functions  are  applied  to  t.he  simulations  [Refs.  36,  57-59). 
a.  Born-Mayer 

The  Born-Mayer  potential  function  is  a  purely  repulsive 
function.  This  function  is  used  for  intermediate  atomic  separations. 
The  form  of  the  function  is 


V(r)  -  a  exp  (-br)  .  (  1  ) 

b.  Mo  1  i  er  e 

The  Moliere  potential  function  is  another  purely  repulsive 
potential  function,  of  a  slightly  different  form.  It  is  called  a 
"Screened  Coulomb  Potential,"  and  is  an  approximation  of  the  Thomas- 
Fermi  screening  function.  The  general  form  of  the  potential  is 
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V ( r )  =  (Z1  Z2  e2  /  r ) [0.35  exp  (-0.3r/a)  +  0.55  exp(-1.2r/a)  + 

0. 1  exp  (-6.0  r/a)  ]  .  (  2  ) 

The  term  "a"  is  called  the  "Firsov  length"  defined  as  follows: 

a  =  K  0.8853  ab  / [Zs  j/2  +  Z2  V2)2/3  .  (  3  ) 

The  "ab"  term  is  the  Bohr  radius, 

ab  z  h2  /  4 tt 2  me2  =  0.5292  A.  (  4  ) 

The  "a"  parameter  may  be  modified,  since  the  Moliere  potential  is  an 
approximation  to  the  Thomas-Fermi  function.  When  K  =  1 ,  the  potential 
is  pure  Moliere;  when  K  is  set  unequal  to  1  then  the  function  is  known 
as  a  "modified  Moliere"  potential  function, 
c.  Morse 

The  Morse  potential  is  both  attractive  and  repulsive, 
depending  upon  the  separation  distance,  "r".  It  has  the  form: 


V (r )  =  De  exp  [-2  a  (r-re)]  -  2De  exp  [-  a(r-re)].  (  5  ) 


The  term,  De  is  the  well  depth,  and  re  is  the  equilibrium  separation  of 
an  atomic  pair,  and  alpha  is  a  scale  factor  which  controls  the  shape  of 
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the  potential  well.  The  function  is  attractive  for  r  >  re,  and  repul- 
s  i  ve  for  r  <  re. 

d.  Composite  Mor se-Mo I i ere 

The  composite  Mor se-Mo 1 i er e  potential  function  is  used  to 
more  closely  model  the  dynamics  of  the  sputtering  events.  The  repulsive 
wall  of  the  Morse  potential  is  joined  with  a  cubic  spline  to  the  Moltere 
potential,  to  form  a  single  potential  function  which  can  be  used  over 
long  ranges.  The  two  functions  are  joined  by  varying  the  alpha  term  of 
the  Morse,  and  the  "a"  term  of  the  Moliere  to  obtain  an  intersection 
which  has  a  smooth,  continuous  slope.  The  repulsive  wall  governs  the 
collision  dynamics,  and  the  attractive  well  of  the  Moliere  controls  the 
sputtering  by  determining  whether  an  atom  will  escape  the  surface  of  the 
target. 

2 .  Selection  of  Potential  Function  Parameters 

a.  Substrate-Substrate,  Adatom-Adatom  Function  Parameters 

A  composite  Mor se-Mo 1 i er e  function  was  used  for  the  solid 
phase  Mo-Mo  and  W-W  potential  functions.  The  specific  parameters  used 
in  this  thesis  follow  those  developed  by  Webb,  with  one  exception.  The 
function  used  for  the  W-W  forces  had  a  discontinuity  as  shown  in  Figure 
5.  This  was  smoothed  by  changing  the  Ra  value.  The  potential  curves 
are  shown  in  Figure  6. 

The  nitrogen-nitrogen  potential  function  is  a  pure  Morse 
potential  function.  The  tabulated  data  for  the  solid  phase  W-W,  Mo-Mo, 
and  N-N  potential  functions  is  listed  below  in  Table  3. 
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o 

o 


Figure  5.  W-W  and  Mo-Mo  Forces 
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POTENTIAL  ENERGY  (EV) 


F i gure  6. 


W-w  and  Mo-Mo 


Potentials  Generated  By  Webb. 
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TABLE  3.  SOLID  PHASE  POTENTIAL  PARAMETERS 


(ev) 

Re 

(A) 

a 

(A"1) 

K 

Ra 

(LU) 

Rb 

(LU) 

Rc 

(LU) 

MO-MO 

0.997 

2.800 

1.519 

0.0 

0.790 

0.830 

2.40 

W-W 

1 .335 

2.894 

1.200 

0.0 

0.80 

1.130 

2.40 

N-N 

7.373 

1.098 

2.700 

0.0 

0.0 

0.0 

1 . 71 

The  a  and  K  terms  were  used  to  match  the  slopes  of  the 
Morse  and  Moliere  potential  functions,  as  mentioned  earlier.  The  value 
of  De  was  selected  in  order  to  obtain  the  proper  cohesive  energy  for  the 
substrate  as  shown  in  Table  1.  The  Re  term  is  the  nearest  neighbor 
distance.  The  spline  boundaries,  Ra  and  R^,  are  in  turn  a  function  of 
the  a  and  K  values  selected,  and  should  be  no  more  than  yt  LU  apart..  This 
topic  will  be  covered  in  more  detail  in  the  following  section.  The 
final  value,  Rc,  is  a  distance  over  which  the  force  and  potential 
calculations  are  made,  a  cut-off  distance.  This  allows  for  interactions 
only  between  nearest  and  next-nearest  neighbors.  This  value  is  set  at 
1.71  LU  for  all  interactions  except  the  substrate-substrate  interac¬ 
tions,  where  it.  is  set.  at  2.40  LU. 

b.  MO-N,  W-N  Potential  Function  Parameters 

Six  new  potential  functions  were  generated  to  examine  the 
effect  of  adatom  placement  on  the  sputtering  cross  section.  The  FORTRAN 
routines  POTTER  and  POTPLOT  were  used  to  help  generate  these  functions. 
The  first  step  was  to  select  an  elevation  to  place  the  adatom,  and  from 
that  to  determine  Re,  the  nearest  neighbor  separation  distance.  Re  is 
the  distance  measured  vertically  between  the  cubic  center  and  the  ada¬ 
tom,  as  shown  in  Figure  7.  These  values  were  then  input  to  POTTER. 
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POTTER  generates  two  data  files  used  by  the  POTLOT  routine 


to  plot  the  potential  and  forces  curves.  The  Ra ,  R^,  k,  and  alpha 
values  are  modified  so  that  a  smooth  curve  is  plotted  for  both  the 
potentials  and  forces.  Once  this  rough  determination  is  completed,  then 
De  is  modified  and  evaluated  by  running  a  single  trajectory  at  one 
time-step  under  CMS  m  order  to  obtain  the  proper  sublimation  energies. 
De  is  modified  so  that  the  sublimation  energy  of  the  nitrogen  to  the 
substrate  is  maintained  at  6.5  eV.  The  binding  energy  of  the  adatom  is 
one-half  the  sublimation  energy,  or  3.25  eV. 

After  the  sublimation  energy  has  been  properly  set,  then 
the  revised  values  for  De  are  once  again  run  through  POTTER  and  POTPLOT 
in  order  to  insure  that  the  curves  have  the  necessary  smooth  slopes. 
These  graphical  methods  for  rapidly  evaluating  potential  function  par¬ 
ameters  have  greatly  eased  the  process  for  generating  potential  func¬ 
tions.  The  values  for  the  Mo-N  and  W-N  potential  functions  used  in  this 
thesis  are  listed  below  m  Tables  4  and  5. 


Figure  7.  Physical  Dimensions  Used  in  Potential  Function 
Calculations. 
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TABLE  4. 


Mo-N  POTENTIAL  FUNCTION  PARAMETERS 


Locat i on 
(A) 

De 

(eV) 

Re 

(A) 

a 

(A)’1 

k 

Ra 

(LU) 

Rb 

(LU) 

Rc 

(LU) 

-0.05 

2.970 

1 .524 

2.60 

0.80 

0.540 

0.560 

1 . 71 

*  # 

0.  146 

2.0428 

1.720 

2.20 

0.80 

0.540 

0.560 

1 . 71 

0.245 

1 . 997 

1 . 738 

2.143 

0.80 

0.500 

0.520 

1 . 71 

*  * 

0.290 

1 . 747 

1 . 864 

2.20 

0.0 

0.560 

0.520 

1 . 71 

0.335 

1 . 672 

1 . 908 

2.20 

0.0 

0.520 

0.560 

1 .71 

0.38 

1 . 60 

1 .950 

2.15 

0.0 

0.720 

0.750 

1 .71 

#  # 

Note 

:  Rows 

marked  with  ##  indicate  parameters 

developed 

by  Webb, 

and  used  m  both 
tion  parameters. 

theses . 

Unmar  ked 

rows  indicate  new 

potential 

func- 

TABLE  5.  W-N  POTENTIAL  FUNCTION  PARAMETERS 


Locat i on 
(A) 

De 

(eV) 

Re 

(A) 

>  Q 

1 

k 

Ra 

(LU) 

Rb 

(LU) 

Rc 

(LU) 

-.05 

3.292 

1.524 

2.85 

0.8 

0.65 

0.67 

1.71  ** 

0. 1623 

2.341 

1 .7448 

2.75 

0.0 

0.75 

0.80 

1 .71 

0.2464 

2.013 

1 .8289 

2.60 

0.0 

0.80 

0.85 

1 . 71 

0.325 

1 . 774 

1 .9075 

2.45 

0.0 

0.85 

0.90 

1 .71 

0.487 

1 .712 

1 .960 

2.32 

0.9 

0.83 

0.85 

1.71  ** 

Note: 

Rows  marked  with  #*  indicate 

parameters  developed  by  Webb, 

and  used  in  both  theses.  Unmarked  rows  indicate  new  potential  func¬ 
tion  parameters. 
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c.  Ion-Substrate,  Ion-Adatom  Functions.  The  potential  func¬ 
tion  parameters  for  the  incident  ion  must  be  considered.  The  target,  is 
bombarded  by  Argon  ions,  but  for  the  purposes  of  these  calculations,  the 
charge  on  the  ion  is  neglected.  The  ion  is  charged  so  that  it  can  be 
accelerated  to  the  desired  energy  to  bombard  the  target.  Since  the 
incident  ion  is  a  noble  gas,  the  event  that  it  will  react  with  the 
molybdenum,  tungsten  or  nitrogen  atoms  other  that  through  simple  colli¬ 
sion  process  is  unlikely.  The  potential  functions  are  correspondingly 
modified  Moliere.  The  potential  function  parameters  are  listed  below  in 
Table  6 . 


TABLE  6.  AR-N,  AR-MO,  AND  AR-W  POTENTIAL  PARAMETERS 


De 

(eV) 

Re 

(A) 

a 

(A)"1 

k 

Ra 

(LU) 

Rb 

(LU) 

Rc 

(LU) 

Ar-N 

0.0 

0.0 

0.0 

0.0 

1  .71 

1  .71 

1  .71 

Ar-Mo 

0.0 

0.0 

0.0 

0.0 

1 .71 

1  .71 

1  .71 

Ar-W 

0.0 

0.0 

0.0 

0.0 

1 . 71 

1  .  71 

1  .71 

Figure  8  shows  the  Nitrogen-Nitrogen  potential  function,  and 
Figure  9  shows  the  Ar-N,  Ar-Mo,  and  Ar-W  potential  functions, 
d.  Vacuum  Phase  Potentials 

Once  atoms  have  sputtered  from  the  surface  of  the  target, 
there  is  a  probability  that  further  interactions  can  occur  in  the  vacuum 
phase.  Vacuum  phase  interactions  were  not  emphasized  in  this  study,  but 
for  the  sake  of  completeness,  the  vacuum  phase  parameters  derived  m 
Webb’s  thesis  are  included  in  Table  7.  There  is  little  available  data 
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Figure  8.  Nitrogen-Nitrogen  Potential  Function. 
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for  the  Mo-Mo,  W-W,  Mo-N,  and  W-N  systems.  The  following  substitutions 


were  made,  using  similar  systems  to  those  examined  in  this 
Ag-Br  system  was  used  fro  Mo-Mo,  and  VO  was  used  for  Mo-N. 
N-N  system  was  available  and  used.  [Refs.  50,  65]. 


study.  The 
Data  on  the 


TABLE  7.  VACUUM  PHASE  POTENTIAL  FUNCTION  PARAMETERS 


De 

(eV) 

a 

(A)"1 

Re 

(A) 

9 

Qx 

Mo-Mo 

1 . 700 

1 . 790 

2.390 

247.000 

0.679 

N-N 

5.  760 

3.359 

1 .098 

2358.570 

14.320 

Mo-N 

6.100 

1 .965 

1 . 524 

1020.000 

4.700 

W-W 

3.  100 

1  .  450 

2.470 

308.399 

0.960 

W-N 

4.900 

2.159 

1 . 738 

967.000 

4.850 

RHC 

2.900 

1 . 560 

2.390 

1050.000 

5.000 

Note:  The 

Rhod i urn 

vacuum  data 

was  used 

i n  s  ome 

cases . 

D.  TARGET  and  IMPACT  AREA  CONCERNS 

One  characteristic  of  Ml  simulations  of  this  type,  is  that  finite 
target  sizes  are  used.  Two  concerns  must  be  balanced  in  determining  the 
optimum  target  size.  First  one  must  try  to  have  a  target  of  sufficient 
size  such  that  the  majority  of  the  sputtering  events  are  contained.  If 
the  target  is  too  small,  then  information  will  be  lost  due  to  a  fa i lure 
of  conta i nment  event.  On  the  other  hand,  an  excessively  large  target 
will  require  excessive  computer  time  to  complete  each  trajectory.  In 
this  thesis,  two  target  sizes  were  used.  At  the  500  eV  energy  level,  a 
19x6x19  target  was  used,  and  at  higher  energy  levels  a  23x8x23  target 
size  was  used.  Webb  used  the  larger  target  in  his  thesis. 

Nitrogen  atoms  were  placed  in  the  four-fold  position  on  both  target 
sizes.  40  atoms  were  placed  on  the  smaller  target,  and  60  atoms  were 


i\lx 


placed  on  the  larger.  In  both  cases  this  resulted  in  a  fractional  cov¬ 


erage  of  approximately  0.4,  which  corresponds  to  the  contracted  domain 
structure  of  Reference  54.  The  nitrogen  adatom  coverage  is  shown  in 
Figure  10,  corresponding  to  the  (001)  face  of  the  substrate  material. 
Table  8  lists  target  area  and  nitrogen  coverage  parameters. 

TABLE  8.  TAR6ET  AREA  AND  ADATOM  COVERAGE  PARAMETERS 


Target  Size 

19x6 

x  19 

23  x  8 

x  23 

Substrate 

Mo 

W 

Mo 

W 

No.  Nitrogen 

40 

40 

60 

60 

Area  (cm^) 

8 . 02x1 0“ 1 4 

8. 1 1x10" 14 

1 . 20x 1 0” 1 3 

1  .21x10" 1 

Tract iona 1 
Coverage 

0.4 

0.4 

0.42 

0.42 

Coverage 

(atoms/cm^) 

4 . 99x 1 0 1 4 

4 . 93x 1 O1 4 

5 . Ox  1 0 1 4 

4 . 96x1 O1 4 

Another  factor  which  impacts  on  the  run  time  of  the  simulation,  and 
the  statistical  accuracy  of  the  results  is  the  number  of  impact  points, 
or  trajectories  run.  The  impact  points  are  designated  by  a  pre-selected 
"pin-pomt",  and  then  the  trajectories  impact  in  a  representative  area 
based  on  this  point.  Two  types  of  pin-points  were  used.  In  one,  the 
"hit"  case,  a  nitrogen  adatom  is  placed  in  the  representative  area,  arid 
can  be  hit  directly  by  the  incoming  ion.  The  other,  "non-hit"  case, 
does  not  have  a  nitrogen  atom  in  the  representative  area.  The  simula¬ 
tions  were  run  in  this  study  with  150  trajectories  per  pin-point,  or  300 
trajectories  per  complete  run.  Previous  studies  used  300  trajectories 
per  pm-pomt,  for  a  total  of  600  trajectories  per  case. 

The  pin-points  for  the  "hit."  case  were  located  at  (8.000,  6.000) 


for  the  19x6x19  target  lattice,  and  (13.000,  15.000)  for  the  23x8x23 


target.  The  non-hit  cases  were  run  with  pin-points  at 


(6.000,  6.000) 


and  (11.000,  11.000) 
x-  and  z-  axes.  The 
i n  F igure  10. 


respect i ve I y . 
r epr esentat ive 


These  pin-pomts  are  measured  on  the 
areas  and  pin-pomts  are  also  shown 


Small  circles  represent  nitrogen  atoms. 

"H"  indicates  "Hit"  representative  area. 
"NH"  indicates  non-hit  representative  area. 


Figure  10.  Top  View  of  Nitrogen  Reacted  Surrace. 
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IV.  RESULTS  AND  DISCUSSION 


A.  GENERAL 

1 .  Spu  t  ter  t  net  or  os  &  Sect  i  ons 

The  concept,  of  the  sputtering  cross  section  is  Key  to  the 
ana  lyses  in  this  study.  The  derivation  of  the  cross  section  is  carried 
out  in  detail  in  Appendix  A.  The  basic  methodology  used  by  Winters  was 
adapted  to  the  simulation,  and  the  following  relation  wus  derived: 

aN  =  -  tin  (1  -  Y/N0) }  A0  .  (  6  ) 

Y  is  the  sputter  yield  from  the  simulation,  is  the  Nitrogen  sputter¬ 
ing  cross  section,  N,-,  is  the  initial  number  of  adatoms  placed  on  the 
target  and  A0  is  the  surface  area  of  the  target. 

2 .  Winters’  Results 

Winters’  experiment  involved  the  bombardment  of  the  (100) 
Plane  of  molybdenum  and  tungsten  with  argon,  xenon  and  helium  ions  in 
energies  from  300  ev  to  5  kev.  The  cases  studied  m  this  thesis  con¬ 
cerned  argon  ions  of  the  500  ev  to  2  kev  range  on  mo (100)  and  W(ioo). 
Winters’  presented  results  from  his  experiments,  and  theoretical  calcu¬ 
lations  from  the  Sigmund-Wi nters  theory  of  multicomponent  sputtering. 
The  results  that  fall  within  the  bounds  of  this  study  are  tabulated 
below  in  Table  9  [Ref.  55  p.  29]. 
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TABLE  9,  WINTERS’  THEORETICAL  AND  EXPERIMENT  CROSS  SECTIONS 
AS  A  FUNCTION  OF  ION  ENERGY 


I  on 

Exper  i merit 

Ca l cu 1  at  ion 

Energy 

(10" 16  cm2) 

( 1 0" 1 6  cm2) 

(ev) 

MO (100)  W(100) 

MO (100)  W(100) 

500 

5.6 

8.6 

6.0  6.9 

1000 

6.5 

8.8 

6,8  T , 4 

2000 

6.8 

9.0 

8.0  8.1 

3000 

6.6 

8.6 

8.9  8.6 

3. 

Previous  Theses 

at  the  Naval 

Postgraduate  School 

a.  Meyer hoff 

Meyerhoff’ s  thesis  in  1983  examined  clean  and  nitrogen 
reacted  molybdenum  and  tungsten  (OOi)  surfaces  bombarded  with  argon  ions 
of  energies  ranging  from  500  ev  to  3  Kev,  Nitrogen  was  placed  in  the 
four-fold  position,  and  was  examined  at  two  elevations  above  the  (001) 
surfaces.  The  effect  of  the  mass  of  the  substrate  was  examined  by 
replacing  the  Mo-mass  with  the  tungsten  mass.  Meyerhoff  concluded  that 
the  sputtering  cross  section  of  the  nitrogen  was  more  dependant  upon  the 
distance  of  the  adatoms  above  the  surface  then  upon  the  mass  of  the  sub¬ 
strate.  [Ref,  49] 

b.  webb 

Webb,  in  early  1966,  continued  the  study  of  nitrogen  on 
molybdenum  and  tungsten.  He  refined  the  potential  functions  used  by 
Meyerhoff,  and  performed  basic  sputtering  research  using  the  qdyn  simu¬ 
lation.  He  studied  the  sputtering  cross  section  of  nitrogen  for  place¬ 
ment  at  three  different  heights  in  relation  to  the  (001)  surfaces  of 
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located  above  the  surface,  as  opposed  to  being  in  the  well  formed  by  the 
s  u  r  f  ace  a  t  onus . 

4 .  Analyses  Conducted  in  This  Study 

The  scope  of  the  analyses  conducted  m  this  study  differed 
somewhat  from  previous  theses.  The  previous  studies  did  not  differen¬ 
tiate  between  the  effects  of  the  incident  ion  initially  hitting  the 
adsorbed  nitrogen,  or  the  substrate.  The  simulations  were  run  with 
impact  points  covering  both  cases,  but  the  results  were  not  looked  at 
separately,  in  this  study,  the  two  cases  were  analyzed  separately. 
This  was  done  m  order  to  gam  insight  into  possible  sputtering  mecha- 
n i sms  as  out  I i ned  ear  1 i er . 

The  nitrogen  sputtering  cross  sections  were  examined  at  500  ev 
(low  energy  sputtering)  and  £  kev  (high  energy  sputtering).  The  hit  and 
no-hit  cases  were  examined  in  order  to  determine  the  sputtering  mecha¬ 
nism  of  the  nitrogen.  The  yield  of  the  substrate  was  examined,  and  the 
ratio  of  the  bare  to  The  reacted  substrate  yields  was  used  as  a  measure 
to  determine  the  relative  effects  of  the  adsorbed  nitrogen.  The  purpose 
of  this  analysis  was  to  determine  if  the  nitrogen  enhanced  the  sputter¬ 
ing  yield  of  the  nitrogen,  particularly  at  the  higher  ton  energies. 

The  effect,  of  the  mass  of  the  substrate  was  examined,  m  order 
to  attempt,  to  resolve  the  difference  between  the  findings  of  Winters  and 
Webb  and  lieyerhoff.  The  technique  used  by  Meyerhoff  was  again  used  m 
this  study,  except  that  the  sputtering  cross  section  of  nitrogen  was 
evaluated  at  three  adatom  placements  as  opposed  to  one,  and  the  hit, 
no-hit  cases  were  examined  separately. 
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B.  MOLYBDENUM  RESULTS 

1 .  Nitrogen  Sputtering  Cross  Sections  and  Adatom  Placement.. 

a.  500  ev  ion  Bombardment 

The  sputtering  of  nitrogen  a?  a  function  of  adatom  place¬ 
ment  was  examined  at  two  energy  levels,  500  ev  and  2  Kev,  in  addition, 
the  "hit"  and  "no-hit"  cases  were  examined  separately.  Table  in  lists 
the  nitrogen  sputter  yield,  and  cross  section  as  a  function  of  adatom 
placement  for  the  500  ev  bombardment.  The  results  are  plotted  in  Figure 
11  for  the  hit,  no-hit  arid  total  cross  sections.  The  plot  of  cross  sec¬ 
tion  vs  atom  placement  results  m  an  apparent,  "well"  corresponding  to  an 
adatom  placement  at  0,245  A  above  the  (100)  surface.  This  was  true  for 
both  the  hit  and  total  cross  section  cases.  The  sputtering  yield  then 
increased  as  a  function  of  adatom  distance  from  the  surface.  The  simula¬ 
tion  cross  section  at  this  point,  is  5,63  x  io-1^  cm2 ,  corresponding  to 
the  experimental  value  obtained  by  Winters  of  5,6  x  io  cm2. 

The  effect,  of  the  ion  hitting  the  nitrogen  is  quite 
apparent.  The  sputtering  cross  section  of  the  hit  case  increased  by  a 
factor  of  from  2  to  4  above  the  non-hit-  case,  and  by  a  factor  of  about. 
1,5  for  the  total  cross  section.  This  indicates  that  the  primary  mecha¬ 
nism  for  the  sputtering  of  the  nitrogen  is  when  the  incident  ion  hits 
the  adatom,  and  it.  m  turn  is  reflected  from  the  substrate,  as  illus¬ 
trated  earlier  in  Figure  2(a).  This  seems  to  be  more  pronounced  as  the 
adatom  is  placed  further  above  the  "well"  at  .245  a.  The  agreement 
between  the  experimental  and  simulation  results  for  this  point  indicates 
that  the  nitrogen  may  be  located  at  approximately  that  height  above  the 
(001)  surface.  The  nitrogen  yields  and  cross  sections  are  1 is  ted  below 
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NITROGEN  SPUTTER  CROSS  SECTIONS 

500  EV  AP.<001>  ON'  MO  (001) 


Figure  11.  Comparison  of  "Hit",  "No-hit."  and  Total  Cross  Sections 
With  Theory  and  Experimental  values  for  500  ev  Ar  on 
Mo l ybdenum. 
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with  H  corresponding  to  the  "hit"  case,  N  for  the  "no-hit"  case 


and  r 


for  the  total  values. 

TABLE  10.  NITR06EN  SPUTTERING  DATA,  Mo(OOI),  500  eV 

Adatom  Sputter  Cross 

Placement  Yield  Section 

(A)  (x  10“ 16  cm2) 


H 

0.41 

9.  1 

0.05 

N 

0.  15 

4.3 

T 

0.28 

6.7 

H 

0.52 

10.5 

0.  146 

N 

0.  1  3 

2.68 

T 

0.33 

6.57 

H 

0.42 

8.47 

0.245 

N 

0.  1  4 

2.81 

T 

0.28 

5.63 

H 

0.63 

12.8 

0.290 

N 

0.20 

4.02 

T 

0.  42 

8.4 

H 

0.66 

13.3 

0.  335 

N 

0.20 

4.02 

T 

8.67 

H 

0.81 

16.3 

0.380 

N 

0.  19 

3.75 

T 

0.50 

10.02 

t>.  2  keV  Bombardment 

Nitrogen  sputtering  when  bombarded  at  2  keV  was  examined 
at  the  same  adatom  locations,  and  the  data  is  tabulated  in  Table  11. 
Figure  12  shows  the  nitrogen  sputtering  cross  sections  as  a  function  of 
adatom  location,  with  the  hit  and  no-hit  case  considered  at  four  loca¬ 
tions.  The  curves  follow  the  same  general  trends  as  Figure  11,  with  some 
exceptions.  The  same  "well"  is  observed,  but  in  this  case,  the  total 
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NITROGEN  SPUTTER  CROSS  SECTIONS 

2.0  KEV  AR<001>  ON  MO(OOl) 


F i gur e  1 2 . 


comparison  of  "H,t",  "No-hit"  ana  Total  cross  sections 
With  Theory  aha  Experimental  values  for  2  k*v  Sr  on 
Mo  I ybdenum. 
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cross  section  of  5.6  x  ICT1*'-'  cnr  is  below  the  experimental  and  theoreti¬ 
cal  values  of  6.8  x  10-1&  cm2  and  8.0  x  10-1&  cm2  respectively.  The 
slope  of  the  sputtering  cross  section  curve  seems  to  level  out  at  adatom 
placements  above  0.29  A.  The  cross  section  for  the  -0.05  placement  is 
below  the  "well”  at  0.245  A.  The  -0.05  point  was  plotted  using  Webb’s 
data.  Both  he  and  heyerhoff  used  14N  in  their  computations,  where  15N 
was  used  in  this  study,  and  in  the  experiments  by  Winters.  The  77.  mass 
difference  may  account  for  part  of  the  difference  in  the  result.  The 
effect  of  the  hit  case  is  again  quite  apparent,  with  significantly 
higher  sputter  cross  sections  resulting  from  the  ion  striking  the 
n i tr ogen  atom. 


TABLE  11.  Mo-N  SPUTTERING  CROSS  SECTIONS,  2.0  KeV 


Adatom 
P I acement 
(A) 


Sputter  Sputter 

Yield  Cross  Section 

(X  10" 16  cm2) 


-0.05 

T 

0.26 

5.2 

H 

0.56 

11.2 

0.  146 

N 

0.15 

3.0 

T 

0.35 

7.01 

H 

0.45 

9.02 

0.245 

N 

0.  1  1 

2.2 

T 

0.28 

5.6 

H 

0.79 

15.9 

0.290 

N 

0.20 

4.0 

T 

0.50 

10.03 

H 

0.80 

16.  1 

0.335 

N 

0.23 

4.60 

T 

0.51 

10.3 

0.380 

T 

0.53 

10.6 
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c.  c omp o r  i  s o n  of  t h*  Sputter  ing  cross  Sections 

The  tot-3 1  sputter  1 09  cross  sections  for  500  ev  000  2  key 
are  plotted  on  the  same  graph  in  Figure  13.  The  general  Trends  are 
apparent.  Both  curves  exhibit,  the  same  "well"  corresponding  to  an  adatom 
placement  of  0.245  a.  The  sputtering  cross  sections  then  increase  as  a 
function  of  increasing  height  of  adatom  placement.  There  seems  to  be 
something  special  about  the  0,245  A  placement.  Since  the  simulation 
cross  section  agreed  so  closely  with  the  experiment,  it  is  a  good  indi¬ 
cation  that  this  might  be  close  to  the  actual  location  of  the  atom.  The 
fact  that  the  2  Kev  curve  so  closely  matches  the  500  ev  curve  gives  fur¬ 
ther  credence  to  this  assumption, 

2 ,  Nitrogen  Sputtering  Cross  Sections  as  a  Function  of  Energy 

The  sputtering  of  nitrogen  as  a  function  of  energy  was 

examined.  The  purpose  of  this  analysis  was  to  compare  the  results  of 

the  simulation  to  the  data  presented  by  Winters  as  listed  in  Table  9. 
Points  were  obtained  at  four  energy  levels,  500  ev  and  1.0,  2,0,  and  3.0 

KeV.  The  analysis  was  conducted  for  adatom  placement  at  0,245  a.  This 

point  was  chosen,  since  it  so  closely  matched  the  experimental  cross 
section  as  mentioned  above.  The  results  are  shown  in  Figure  1 4 ,  with 
the  corresponding  experimental  and  theoretical  values  from  Winters’ 
paper,  The  values  at  500  ev  and  3.0  kev  are  within  the  bounds  of  the 

experimental  and  theoretical  values,  but  the  points  at  1  Kev  and  2  Kev 

are  beyond  the  limits.  The  variation  of  the  simulation  results  from 
those  of  Winters  is  a  measure  of  the  uncertainty  of  the  simulation 

model .  However,  the  general  trend  of  the  simulation  cross  sections  d<> 

correlate  with  Winters’  results. 
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Figure  13.  Comparison  of  Total  Cross  Sections  for  500  ev  and  2  keV  Ar 
on  Molybdenum. 
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NITROGEN  SPUTTER  CROSS  SECTIONS 

N  AT  0.245  AN.  MO(OOi) 


Figure  14.  Nitrogen  sputtering  cross  Section  as  a  Function  of  Energy 
for  Adatoms  Placed  at  0.24?  a  on  Mo(ioo), 
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3, 


Spu 1 1 er  « fVj  of  The  Substrate 


The  ratio  of  the  reacted  substrate  yield  to  the  bare  substrate 

« 

yield  was  examined  at  500  ev  and  2  kev.  The  hit,  no-hit.  and  total 
yields  were  considered.  The  results  are  listed  m  Table  12,  arid  are 
plotted  in  Figure  15.  in  all  cases,  for  ail  adatom  placements  examined, 
the  yield  of  the  reacted  surface  was  lower  than  that  of  the  bare  sur¬ 
face,  Add i t i ona 1 1 y ,  the  substrate  yields  of  the  hit  case  were  lower 
than  for  the  total  and  non-hit  case.  This  indicates  that  the  nitrogen 
adatoms  absorb  some  of  the  energy  of  the  incident  particle,  and  there¬ 
fore  less  momentum  is  available  to  initiate  the  collision  cascade  in  the 
substrate. 

The  case  was  the  same  at  2  kev.  The  sputtering  yields  of  the 
substrate  were  lower  in  each  case  when  the  adatom  was  initially  struck 
by  the  ion.  The  incident  ion  at  2  kev  still  failed  to  have  sufficient 
energy  to  impart  enough  momentum  to  the  adatom  for  it  to  act  as  a  second 
projectile  and  thereby  enhance  the  sputter  yield  of  the  substrate.  The 
yields  of  the  substrate  and  the  ratios  with  the  bare  substrate  for  2  kev 
are  also  listed  in  Table  12. 

C,  TUNGSTEN  RESULTS 

1 .  Nitrogen  sputtering  cross  Sections  and  Adatom  Placement 

a.  500  ev  ion  Bombardment 


The 

sputter i ng 

of  nitrogen  from  (001)  tungsten 

was 

examined 

at  500  ev 

and  2  kev. 

The  hit  and 

no-hit  scenarios 

were 

exami ned 

as  with 

mo l ybdenum. 

The  results  are 

listed  in  Table  13, 

and 
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TABLE  12. 


SPUTTER  YIELD  OF  SUBSTRATE,  Mo  (100.) 


500 


kev 


Adatom 


P I  a cement 

(A) 

Yield 

Wyb 

Yield 

Wyb 

H 

3.16 

NA 

Bare 

N 

3.41 

NA 

T 

3.28 

NA 

3,02 

NA 

H 

1 . 99 

0, 630 

-.05 

N 

2.29 

0,671 

T 

2,14 

0 , 652 

2.29 

0,76 

H 

1 .81 

0,573 

2.01 

0.67 

0.146 

N 

2.37 

0.695 

3.07 

1.02  # 

T 

2.09 

0 . 637 

2,54 

0,84 

H 

1 .91 

0, 604 

1  .  65 

0.61 

0.  245 

N 

2.45 

0.718 

2.66 

0,88 

T 

2,  17 

0,661 

2.25 

0.74 

H 

2 . 03 

0,642 

2.28 

0.75 

0.290 

N 

2.51 

0 , 736 

3.15 

1.04  * 

T 

2.27 

0.692 

2,72 

0,90 

H 

2.08 

0,  658 

2,09 

0, 69 

0,335 

N 

2.60 

0.762 

3.  17 

1,05  * 

T 

2,34 

0.713 

2 , 63 

0,87 

H 

2.03 

0.642 

0.380 

N 

2.59 

0 . 760 

T 

2,31 

0,704 

2,25 

0,74 

*  *  It 

shou  l  d 

be  noted  that 

the  results 

for  the  bare 

substrate 

Mo- bare  case  was  not 
for  the  ratio  of 


were  taken  from  reference  50.  The  2  kev 
re-run.  This  could  be  part  of  the  reason 
sputter  yields  for  the  new  potential  functions  to  be  greater 
than  1.  it  does  not  seem  reasonable  that,  the  yield  of  the 
reacted  substrate  should  be  greater  than  the  bare  substrate 
when  the  adatom  is  not  struck. 
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RATIO  OF  REACTED  TO  DARE  YIELDS 

0.2  0.4  0.0 


SPUTTER  YIELDS  OF  SUBSTRATE 

P.EACTED/BARE  SPUTTER  YIELD,  MO 


Figure  15. 


Ratio  of  the  Sputter  Yield  of  the  Reacted  Substrate  to  the 
Yield  of  the  Bare  Substrate. 
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are  p  lotted  in  Figure  if’-.  The  same  genera I  pattern  wa^  Fyur.  <3  f<,.r 
molybdenum.  The  sputtering  cross  section  for  the  hit  scenar 10  is  ago  in 
significantly  greater  than  for  the  no-hit  case,  for  all  adatom  place¬ 
ments.  Therefore,  as  with  molybdenum,  the  majority  of  the  nitrogen  is 
sputtered  as  a  result  of  a  direct  collision  by  the  incident  ion.  The 
total  sputtering  cross  sections  generally  follow  the  experimental  val¬ 
ues  until  the  point  corresponding  to  the  0.2464  a  adatom  location,  where 
the  slope  of  the  cross  section  plot  rapidly  increases. 

b .  2  k  e v  I o n  Bomba  r dme n t 

The  sputtering  cross  section  and  yield  data  for  the  2  kev 
bombardment  is  sumnar i zed  in  Table  m  and  plotted  in  Figure  it.  There 
is  no  "well"  as  was  observed  for  the  molybdenum  case,  however  there  is  a 
break  m  the  slope  of  the  curve  at  the  point  corresponding  to  adatom 
placement  at  0.2464  a.  While  not  plotted,  the  sputtering  cross  sections 
for  the  hit  case  were  again  significantly  greater  in  all  instances  than 
for  the  non-hit  case. 

c .  C ornp a r  i  s o n  of  S p u 1 1 e r  i  n g  C r o s s  S e c 1 1  o n s 

The  total  nitrogen  sputtering  cross  sections  for  the  500 
ev  and  2  kev  bombardments  are  plotted  in  Figure  19.  The  two  curves  fol¬ 
low  the  same  general  trends.  A  "well"  is  not  observed,  but  there  is  an 
inflection  point  in  the  curves  at  a  point  corresponding  to  an  adatom 
Placement  of  0.2464  A.  One  special  note,  the  sputtering  cross  section 
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CROSS  SECTION  (X  E-16  SQ  CM) 

6.0  0.0  10.0  1L\0  11.0  16.0  10.0  20.0  22.0 


F i gur  e 


NITROGEN  SPUTTER  CROSS  SECTIONS 

500  EV  AP.<001>  OX  ff(OOl) 


16.  Comparison  of  Hit,  No-hit  ana  Total  Nitrogen  Sputtering 

cross  Section  with  Experiment  aria  Theory  for  too  ev  Argon 
on  W(ioo)  ' 
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NITROGEN  SPUTTER  CROSS 'SECTIONS 


2  KEV  AR<001>  ON  ff(OOl) 


Figure  17.  Nitrogen  sputtering  Cross.  Section  as  a  Function  of  Adatom 
Placement  for  2  kev  Argon  on  W(IOO). 
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CROSS  SUCTION  (X  E-1G  SQ  CM) 


NITROGEN  SPUTTER  CROSS  SECTIONS 

.5  AND  2  KEV  AR<001>  ON  1V(001) 


Figure  18.  Comparison  of  Total  Nitrogen  Sputtering  Cross 
Sections  for  500  eV  and  2  keV  Argon  on  W(100). 
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TABLE  13.  WtlOO)  NITROGEN  SPUTTERING  RESULTS,  500  ev' 


Adatom  Sputter  Sputter 


P 1 acement 
(A) 


H 

-.05  N 

T 

H 

0.1623  N 

T 

H 

0.2464  N 

T 

H 

0,325  N 

T 

H 

0.467  N 

T 


Y  i  e  l  d  c  r  o  5  *  s  e  c  t  i  o  n 

(x  10" 16  cm2) 


0.60 

12,3 

0,  23 

4.68 

0,42 

l r> 

CO 

0, 69 

14,1 

0.23 

4 . 68 

0.46 

9 , 36 

0. 68 

13,9 

0.27 

5,49 

0.47 

9 , 59 

0,63 

17.0 

0. 32 

6,51 

0,57 

11,6 

0, 99 

20,3 

0.33 

6,72 

0.66 

13,5 

TABLE  14.  NITROGEN  YIELD  AND  SPUTTERING  CROSS  SECTIONS  FROM  TUNGSTEN 

2.0  kev 


Adatom  Total 

Placement  Sputter 

(A)  Yield 


Sputter 
Cross  Section 
(x  10~1Fj  cm2) 


.05 

0.28 

5,  67 

,  1623 

0,35 

7.08 

,  2464 

0,36 

7.28 

,325 

0.40 

6,09 

.  487 

0,  45 

9,  1  1 
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Of  nitrogen  at  500  eV  i  *  s  i  90  1  f  1  Cant  ly  0  f  e  a  ter  fhan  at  d  KeV  for  jl  I 
add  tom  placements .  This  is  similar  to  the  findings  in  the  two  previous 
theses  [Ref.  49,  50],  This  conflicts  with  the  experimental  and  theoret¬ 
ical  values  presented  t>y  Winters.  This  is  an  indication  that  the  model 
used  m  the  simulation  is  not  correctly  portraying  the  behavior  of  the 
n 1 trogen-tungsten  system. 

£,  Nitrogen  sputtering  as  a  Function  of  Energy 

In  an  effort  to  find  an  adatom  placement  location  that  agreed 
with  the  experimental  results,  the  total  nitrogen  sputtering  cross  sec¬ 
tion  was  plotted  as  a  function  of  energy  for  all  five  adatom  locations 
tested.  The  results  are  shown  in  Figure  19,  In  all  cases,  there  was  a 
strong  negative  slope,  that  did  not  correspond  well  with  the  slight  pos¬ 
itive  slope  of  the  experimental  and  theoretical  values.  Again  this  cor¬ 
responds  to  Webb’s  and  Meyer hoff’s  findings.  This  indicates  that  the 
potential  function  may  not  be  reflecting  the  true  behavior  of  the 
nitrogen  on  the  tungsten. 

3.  Sputtering  of  the  Substrate 

The  analysis  of  the  ratio  of  the  sputtering  yield  of  the 
reacted  substrate  to  the  bare  substrate  was  not  conducted  to  the  same 
detail  as  with  molybdenum.  The  results  for  the  cases  examined  are  sum¬ 
marized  below  in  Table  15,  in  all  cases,  with  the  exception  of  the  hit- 
case  for  placement  below  the  surface  at  500  ev,  the  yield  of  the  reacted 
substrate  was  lower  than  that  of  the  bare  substrate.  This  wa*  the  case 
both  at  500  ev  and  2  KeV.  This  indicates  that  for  these  two  levels,  the 
adatoms  again  reduce  the  energy  of  the  ion  such  that  less  momentum  1  s 
transferred  to  the  substrate,  and  the  yield  of  the  substrate  is  lowered. 
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NITROGEN  SPUTTER  CROSS  SECTIONS 


Figure  19.  Comparison  of  Nitrogen  Sputtering  Cross  sections  as 
a  Function  of  Energy  and  Adatom  Location,  W(100). 
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One  interesting  observation  i  that  the  yield  of  the  substrate  a  t  2  Key 
did  i nor  ease,  a s  expected,  at  the  higher  energy.  This  despite  the  fact 
that  the  simulation  cross  section  of  nitrogen  sputtering  was  higher  at 
500  ev  than  at  2  Kev. 

TABLE  15.  SPUTTER  YIELD  OF  SUBSTRATE  W(100) 

500  ev  2  Kev 

Adatom 

Placement  Yield  yn/yb  Yield  yn/yb 

A 


H 

2.67 

Bare 

N 

3.26 

NA 

T 

c ,  96 

2 . 63 

NA 

H 

1 . 79 

0. 67 

2.16 

0.82 

-0.05 

N 

1.75 

0,53 

1 ,57 

0,  60 

T 

1 . 77 

0. 59 

1 .87 

0.71 

H 

1 .83 

0.68 

0. 1623 

N 

1 .65 

0.57 

T 

1.85 

0.62 

2.06 

0.78 

H 

1 . 77 

0. 67 

0.2464 

N 

1 ,9 

0.57 

T 

1 .84 

0.62 

2.35 

0.89 

H 

1 . 76 

0. 66 

0.325 

N 

2.05 

0.63 

T 

1 .90 

0.64 

2.50 

0.95 

H 

1.73 

0.65 

0.467 

N 

2.23 

0.68 

T 

1.98 

0 . 66 

1 .96 

0.74 

COMPAR 1  SON 

OF 

MOLYBDENUM 

AND  TUNGSTEN 

RESULTS 

1 .  Compar 

ison  of  Cross 

Sect i ons 

The  sputtering  cross  sections  of  nitrogen  from  mol y Oden urn  and 


tungsten  are  plotted  m  Figure  20,  with  the  adatom  placement  normalized 
to  lattice  units.  The  molybdenum  curves  exhibit  the  "well"  at.  0.156  lu, 
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CROSS  SECTION  (X  E-1G  SQ  CM) 

3.0  4.0  5.0  G.O  7.0  0.0  9.0  10.0  11.0  1Z 


NITROGEN  SPUTTER  CROSS  SECTIONS 

.5  AND  2KEV  AR  ON  W  AND  MO(OOl) 


Figure  20.  Comparison  of  Nitrogen  Sputtering  Cross  Sections 

for  Mo(IOO)  and  W(100)  at  500  eV  and  2  keV,  Adatom 
Placement  Normalized  to  Lattice  Units. 
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but  the  tuny* ten  curves  show  a  definite  "break"  in  slope  at  the  . . 
relative  point.  This  point  corresponds  to  0.245  a  for  molybdenum,  and 
0.2464  a  for  tungsten.  The  potential  function  for  the  0.2464  LU  place¬ 
ment  was  derived  so  that  the  behavior  at  the  point  correspond i ng  to  the 
molybdenum  well  could  be  examined.  Again  the  major  difference  between 
the  two  substrates  is  the  fact  that,  the  sputtering  cross  section  at  500 
ev  is  higher  than  at  2  kev  for  tungsten. 

The  hit  and  no-hit-  analysis  was  conducted  for  tungsten  also, 
and  the  same  general  trend  was  noted.  Significantly  higher  sputtering 
cross  sections  were  noted  for  the  hit  scenario.  This  was  true  at  all 
energy  levels  and  adatom  placements  examined.  This  indicates  that  the 
primary  mechanism  for  the  sputtering  of  nitrogen  is  due  to  direct  colli¬ 
sion  by  the  incident  ion, 

2 .  Comparison  of  Substrate  Sputtering  Yields 

The  ratios  of  the  sputtering  yields  of  the  reacted  substrates 
to  the  bare  substrates  at  500  ev  and  2  kev  bombardments  indicate  that 
the  nitrogen  does  not  enhance  the  sputtering  of  the  substrate.  This  is 
observed  even  at  the  relatively  high  energy  level  of  2  kev.  Both 
mo l ybdenum  and  tungsten  exhibited  this  behavior.  The  incident  ion,  even 
at  2  kev  failed  to  transfer  enough  momentum  to  the  nitrogen  for  it  to 
act  as  a  second  projectile  to  enhance  the  sputtering  of  the  substrate. 

E.  DETERMINATION  OF  POSSIBLE  MASS  EFFECTS 

One  of  the  major  conclusions  in  Winters’  earlier  works  was  the  fact 
that  the  sputtering  cross  section  seemed  to  be  greater  for  the  substrate 
with  the  higher  mass  [Ref,  49,  55],  Meyerhoff  and  Webb  examined  this 
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j  5  p  r  r  t.  ,  a  0  '"I  C  0  0  C  I  U  0  0  ('1  t  h  6 1  a  d  3  1 010  p  I  ,j  0  er fir  r  1 1  had  d  ‘J  r  e  i  f  « r  ■  £  f  f  e  f  .  T  h  c 

mass  effect  was  examined  in  this  thesis  by  selecting  three  adatom  loca¬ 
tions  m  the  vicinity  of  the  "well"  at-  0.245  a.  for  the  Mo (100)  surface, 
and  replacing  the  mass  of  the  molybdenum  with  the  tungsten  mass,  leaving 
all  other  aspects  of  the  potential  functions  the  same.  The  results  ore 
plotted  in  Figure  21,  for  the  hit,  no-hit  and  total  cross  section  cases, 
in  all  cases,  the  sputtering  cross  sections  were  higher  for  the  higher 
mass.  This  indicates  that  the  mass  of  the  substrate  does  have  on  effect 
on  the  sputtering  cross  section  of  the  nitrogen. 

Figure  22  shows  the  total  nitrogen  sputtering  cross  sections  for 
500  ev  mo(IOO),  Mo  * ( 1 00)  (the  high-mass  molybdenum),  and  w(ioo).  The 
graph  is  normalized  for  adatom  placements  at  approximately  corresponding 
locations.  The  trend  for  higher  sputter  cross  sections  for  the  higher 
mass  is  obvious.  The  heavy  Mo*  has  higher  nitrogen  sputtering  cross  sec¬ 
tions  than  the  true  molybdenum,  and  the  tungsten  has  the  highest  cross 
sections  of  the  three.  Through  the  use  of  the  computer  simulation,  the 
mass  effects  have  been  able  to  be  truly  isolated,  more  so  than  is  pos¬ 
sible  in  a  physical  system.  This  is  one  valuable  aspect  of  simulations. 

One  possible  reason  for  the  higher  sputtering  cross  section  for  the 
heavier  substrate  is  that  the  greater  mass  of  the  substrate  provides  a 
better  "springboard"  for  the  nitrogen  atoms.  When  the  nitrogen  atoms 
are  Knocked  from  their  positions  by  the  incident  ions,  they  relatively 
light  nitrogen  atoms  will  rebound  with  more  initial  momentum  conserved 
from  the  heavier  tungsten  atoms  than  from  the  lighter  molybdenum  atoms. 

This  model  for  the  increase  in  sputtering  yield  with  higher  sub¬ 
strate  mass  follows  if  the  primary  mechanism  for  the  sputtering  of  the 
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CROSS  SECTION  (X  E-1G  SQ  CM) 

2.0  4.0  0.0  0.0  10. 0  12.0  14.0  16.0  10.0  20.0 


NITROGEN  SPUTTER  CROSS  SECTIONS 

500  EV  AR  ON  MO  AND  MO "(001) 


Figure  21.  Comparison  of  Nitrogen  Sputtering  Yields  for 
Mo ( 1 00)  and  Mo#(100),  500  eV  Argon. 
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CROSS  SECTION  (X  E-1G  SQ  CM) 


NITROGEN  SPUTTER  CROSS  SECTIONS 


ADATOM  PLACEMENT  (LU) 


Figure  22.  Comparison  of  Total  Nitrogen  Sputtering  Cross 

Sections  for  Mo(IOO),  Mo* (100),  andW(IOO),  500  eV 
Argon. 
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nit  rogen  is  d  i  r  ec  t  COM  i  s  i  o  o  Oy  the  incident  ion.  if  the  n  1 1  r  o  g  e  n  atoms 
wer e  knocked  off  by  reflected  ergon  ions,  or  by  sputtering  substrate 
material s  to  a  high  degree,  the  effect  of  the  mass  of  the  substrate 
would  seem  to  be  less  important. 

F.  COMPARISON  WITH  PREVIOUS  SIMULATIONS 

The  results  of  this  study  indicates  that  there  is  a  significant 
mass  effect,  this  is  contrary  to  the  findings  of  Webb  and  Meyerhoff. 
The  earlier  studies  did  not  isolate  the  effects  of  the  direct  interac¬ 
tion  of  the  ion  with  the  adatom  to  the  degree  that  it.  was  done  in  this 
study.  The  trend  for  the  higher  nitrogen  sputtering  cross  section  at 
the  lower  energy  for  tungsten  is  consistent  through  the  three  studies. 

The  other  aspects  of  the  computer  simulation,  such  as  the  angular 
and  energy  distribution  plots  remained  consistent  with  earlier  works. 
For  a  me* re  detailed  analysis  of  the  "traditional"  sputtering  simulation 
analysis,  the  reader  is  referred  to  References  49  and  50. 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  results  of  this  simulation  confirm  the  findings  of  Webb  and 
Meyerhoff  that  the  placement  of  the  adatoms  is  most  likely  above  the 
crystal  surface.  A  location  of  approximately  0.245  A  for  the  molybdenum 
surface  results  in  a  nitrogen  sputtering  cross  section  that  agrees  with 
the  experimental  results  of  Winters. 

Comparison  of  the  "hit"  and  "no-hit"  scenarios  indicates  that  the 
primary  mechanism  for  the  sputtering  of  the  adatom  is  due  to  direct 
collisions  with  the  incident  ion.  In  all  cases  examined,  the  sputtering 
cross  sections  derived  when  the  adatom  was  struck  by  the  incident  ion 
were  larger  than  in  the  cases  when  they  were  not  struck.  This  was  found 
at  both  energy  levels  examined,  500  eV  and  2  keV;  and  for  both  Mo (100) 
and  W ( 1 00) .  This  agrees  partially  with  Winter’s  conclusions,  but  does 
not  show  a  significant  increase  in  non-hit  sputtering  at  higher  energies 
he  proposed. 

The  mass  effect  of  the  substrate  was  examined,  and  it  was  deter¬ 
mined  that  the  sputtering  cross  section  of  the  adatom  was  enhanced  with 
the  increase  of  mass  of  the  substrate,  when  all  other  factors  are  held 
constant.  This  contradicts  Meyerhoff’s  findings,  but  only  one  point  was 
considered  in  his  analysis. 

The  ratio  of  the  sputtering  yield  of  the  bare  substrate  to  the 
reacted  substrate  was  evaluated  for  both  systems.  The  analyses  indi¬ 
cated  that  the  adatoms  decrease  the  amount  of  energy  imparted  to  the 
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substrate,  and  reduces  the  sputter  yield  for  the  "hit" 


case. 


Th  i  s 


momentum  loss  was  observed  at  500  eV  and  2  KeV.  This  process  cocid 

change  at  higher  energy  levels  (greater  than  2  keV)  of  the  incident  par¬ 
ticles. 

A  significant  difference  was  found  for  tungsten  between  the  simula¬ 
tion  and  experimental  and  theoretical  sputtering  cross  sections.  The 
nitrogen  sputtering  cross  section  for  the  500  eV  bombarding  energy  was 
higher  than  for  the  2  kev  energy,  which  was  in  contradiction  to  Winter’s 
results.  A  similar  finding  was  made  by  Webb.  This  indicates  that  per¬ 
haps  the  pair  potentials  used  to  model  the  nitrogen  on  tungsten  are  not 
correct.  It  has  been  determined  that  the  oxygen-tungsten  system  behaves 
quite  differently  than  the  oxygen-molybdenum  system,  so  perhaps  this 
trend  continues  for  the  Mo-N  and  W-N  systems  [Ref.  67]. 

Further  simulation  studies  should  be  made  to  complete  the  compari¬ 
son  with  Winters’  study.  Studies  should  be  made  using  helium  and  xenon 
ions  in  order  to  investigate  the  effect  of  the  mass  of  the  incident  ion. 
Additional  studies  should  be  conducted  on  tungsten  to  determine  if  the 
lower  sputtering  cross  section  at  2  keV  was  due  to  nitrogen  atoms  being 
captured  as  interstitials  in  the  tungsten  lattice,  or  if  there  is  a 
defect  in  the  potential  functions.  Higher  energy  studies  should  be 
conducted  to  identify  the  point  at  which  the  sputtering  yield  of  the 
substrate  will  be  enhanced  by  the  adatoms. 
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APPEND  I  X 


DERIVATION  OF  SPUTTERING  CROSS  SECTIONS 

The  methodology  used  for  the  derivation  of  the  sputtering  cross 
section  in  this  thesis  is  based  on  the  technique  used  by  Winters  in 
reference  55.  The  notation  used  in  this  appendix  closely  parallels  that 
of  W i nter s . 

The  rate  that  nitrogen  is  sputtered  from  the  surface  of  the  target 
is  given  by  the  relation: 

d©N 

-rn  -  -  -  -  a  eN  V+  t  ,  (  7  ) 

dt 


where  the  following  are  defined: 

Rn  =  the  nitrogen  sputtering  rate,  {atoms  /  cm2  -  secj; 

©N  =  the  nitrogen  atom  concentration  at  time  =  t,  {atoms  /  cm2]; 

a  =  the  sputtering  cross  section,  {cm2]; 

V+  =  the  incident  ion  flux,  {ion  /  cm2  -  sec]; 
t.  =  the  time  in  seconds,  {sec};  and 
V+t  =  the  ion  dose,  (fluence),  {ion  /  cm2}. 


Solving  equation  7  yields, 

In  (©N  /  ©No)  =  -ay+ t,  (  8  ) 

where  ©n0  =  the  initial  nitrogen  concentration  at  t  =  0.  Rearranging 
equation  8,  and  solving  for  a  results  in  the  following  relation  for  the 
cross  section; 

a  ^  -in  (©N/eNo)  /  V  +  t  .  (  9  ) 
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The  ©n/©n0  t:erm  ,s  next  examined. 


The  following  relations  are 


def i ned, 


No 

NTRAJ 

NTRAJ 

NSPUTT 

A0 


the  number  of  nitrogen  atoms  on  the  surface  at  t  =  0, 

the  number  of  trajectories  run,  also 

the  number  of  incident  ions  on  the  target, 

the  number  of  nitrogen  atoms  sputtered  from  the  surface, 

the  area  of  the  target  {cm2). 


Consequently,  for  each  trajectory  we  have 


0N  =  (N0  -  NSPUTT)  /  A0, 


and  ©no  -  N0/A0- 

Thus , 

eN 

(N0  - 

NSPUTT)  /  A0 

(N0  -  NSPUTT) 

eNo 

^0  /  ^0 

No 

Surrming  over  all  trajectories, 


©N  ^NTRAJ 

©No  '  =  1 


(N0  -  NSPUTT)  j 

No 


=  1  - 


NTRAJ 


NSPUTT, 


(N0)  (NTRAJ) 


NTRAJ 

If  we  define  £  NSPUTT  -  TOTSPUT,  then 

i  =  1 
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1 


TOTSPUT 


(  10  ) 


©N 

®No 


(N0)  (NTRAJ) 


Incorporating  the  relation  for  the  beam  fluence,  from  equation  9 
yields  the  following  relation  for  the  cross  section; 


but 


o  = 


TOTSPUTT 


NTRAJ 


TOTSPUTT 


(N0)  (NTRAJ) 


) 


=  Y I  ELD,  so  as  a  r  esu  1 1; 


A 


O' 


Equation  11  is  the  relation  used  to  calculate  the  sputtering  cross 
section  for  the  adsorbed  nitrogen  from  the  tungsten  and  molybdenum  sur¬ 
faces  . 
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